clear
set more off

*ssc install blindschemes /* requires installing plotplain scheme - type findit plotplain in command window*/

*** merge SSI education data ***

insheet using "ssi-education.csv", comma name
rename IdParameter psid
save "ssi-education.dta", replace

*** STUDY 1, WAVE 2 ***

use "Zika_CT_experimentwave_2.dta", clear

/* Recode Wave 2 variables */
recode Q340 (10=1 "Not accurate") (11=2 "Somewhat accurate") (12=3 "More or less accurate") (13=4 "Very accurate") (-99=.), gen(zika_sex_w2)
/* A person can become infected with Zika by having sex with someone who is
infected. */
/*
8		not recorded			-99
439		Not accurate			10
128		Somewhat accurate		11
188		More or less accurate	12
193		Very accurate			13
*/

recode Q342 (10=4 "Not accurate") (11=3 "Somewhat accurate") (12=2 "More or less accurate") (13=1 "Very accurate") (-99=.), gen(zika_handshake_w2)
/* A person can become infected with Zika through casual contact with someone 
who is infected. */
/*
5		not recorded			-99
819		Not accurate			10
63		Somewhat accurate		11
45		More or less accurate	12
24		Very accurate			13
*/

recode Q344 (6=1 "Not accurate") (10=2 "Somewhat accurate") (11=3 "More or less accurate") (12=4 "Very accurate") (-99=.), gen(zika_bite_w2)
/* A person may be infected with Zika by the bite of a mosquito carrying the
virus. */
/*
5		not recorded			-99
14		Not accurate			6
29		Somewhat accurate		10
45		More or less accurate	11
863		Very accurate			12
*/

recode Q346 (10=1 "Not accurate") (11=2 "Somewhat accurate") (12=3 "More or less accurate") (13=4 "Very accurate") (-99=.), gen(zika_gmo_w2)
/* Genetically modified mosquitoes - not wild mosquitoes - caused Zika's
 outbreak in Latin America. */
/* 
3		not recorded			-99
465		Not accurate			10
209		Somewhat accurate		11
186		More or less accurate	12
93		Very accurate			13
*/

recode Q348 (10=1 "Not accurate") (11=2 "Somewhat accurate") (12=3 "More or less accurate") (13=4 "Very accurate") (-99=.), gen(zika_neuro_w2)
/* Being infected with the Zika virus increases the likelihood of developing 
neurological problems. */
/*
7		not recorded			-99
87		Not accurate			10
124		Somewhat accurate		11
292		More or less accurate	12
446		Very accurate			13
*/

recode Q354 (10=1 "Not accurate") (11=2 "Somewhat accurate") (12=3 "More or less accurate") (13=4 "Very accurate") (-99=.), gen(micro_larvicide_w2)
/* The prevalence of microcephaly in Brazil has increased as a result of 
larvicides added to water to kill mosquitoes */
/*
3		not recorded			-99
589		Not accurate			10
167		Somewhat accurate		11
111		More or less accurate	12
73		Very accurate			13
*/

recode Q356 (10=1 "Not accurate") (11=2 "Somewhat accurate") (12=3 "More or less accurate") (13=4 "Very accurate") (-99=.), gen(micro_vaccine_w2)
/* The prevalence of microcephaly in Brazil has increased as a result of the
government determining increased use of vaccines among pregnant women. */
/*
3		not recorded			-99
625		Not accurate			10
152		Somewhat accurate		11
92		More or less accurate	12
71		Very accurate			13
*/

recode Q364 (4=1) (5=2) (6=3) (7=4) (8=5) (9=6) (10=7) (11=8) (12=9) (13=10) (-99=.), gen(policy_dtap_w2)
/* Do you approve or disapprove of the government's recommendation that pregnant
women receive the DTaP vaccine? */
/* 1 = strongly disapprove, 10 = strongly approve */
/*
1	not recorded		-99
24	strongly disapprove	1
13						2
13						3
28						4
98						5
55						6
93						7
164						8
92						9
346	strongly approve	10
*/

recode Q362 (4=1) (5=2) (6=3) (7=4) (8=5) (9=6) (10=7) (11=8) (12=9) (13=10) (-99=.), gen(policy_larvicide_w2)
/* Do you approve or disapprove of the use of larvicides in water to prevent 
mosquito breeding in areas where Zika is believed to be present? */
/* 1 = strongly disapprove, 10 = strongly approve */
/*
2	not recorded		-99
21	strongly disapprove	1
12						2
17						3
33						4
79						5
58						6
96						7
154						8
94						9
361	strongly approve	10
*/

recode Q366 (4=1) (5=2) (6=3) (7=4) (8=5) (9=6) (10=7) (11=8) (12=9) (13=10) (-99=.), gen(policy_gmo_w2)
/* Do you approve or disapprove of the release of genetically modified 
mosquitoes to the environment to combat the spread of Zika? */
/* 1 = strongly disapprove, 10 = strongly approve */
/*
0	not recorded		-99
41	strongly disapprove	1
11						2
22						3
29						4
90						5
64						6
95						7
165						8
86						9
324	strongly approve	10
*/

recode Q368 (4=1) (5=2) (6=3) (7=4) (8=5) (9=6) (10=7) (11=8) (12=9) (13=10) (-99=.), gen(policy_homesearch_w2)
/* Do you approve or disapprove of government health workers entering private 
homes to look for mosquito breeding points / outbreaks when the resident or 
owner is not present? */
/* 1 = strongly disapprove, 10 = strongly approve */
/*
2	not recorded		-99
71	strongly disapprove	1
17						2
21						3
23						4
80						5
63						6
80						7
130						8
84						9
356	strongly approve	10
*/

recode Q372 (7=1 "No risk") (8=2 "Very low risk") (9=3 "Low risk") (10=4 "Moderate risk") (12=5 "High risk") (13=6 "Very high risk") (-99=.), gen(risk_human_w2)
/* How much risk do you think the Zika virus presents to human health in
Brazil? */
/*
2	not recorded	-99
2	no risk			7
16	very low risk	8
32	low risk		9
191	moderate risk	10
351	high risk		12
332	very high risk	13
*/

recode Q380 (7=1 "No risk") (8=2 "Very low risk") (9=3 "Low risk") (10=4 "Moderate risk") (12=5 "High risk") (13=6 "Very high risk") (-99=.), gen(risk_personal_w2)
/* How much risk do you believe the Zika virus represents to you personally? */
/*
2	not recorded	-99
31	no risk			7
87	very low risk	8
168	low risk		9
312	moderate risk	10
198	high risk		12
128	very high risk	13
*/

recode Q384 (1=1 "Not at all concerned") (2=2 "Not very concerned") (3=3 "Somewhat concerned") (4=4 "Very concerned") (5=5 "Extremely concerned") (-99=.), gen(concern_w2)
/* How concerned are you personally about the potential health risk for Zika? */
/* 
11	not recorded			-99
29	not at all concerned	1
110	not very concerned		2
303	somewhat concerned		3
281	very concerned			4
192	extremely concerned		5
*/

recode Q390 (4=5 "Always") (5=4 "Almost always") (6=3 "Sometimes") (7=2 "Rarely") (8=1 "Never") (-99=.), gen(protect_sleeves_w2)
/* How often do you protect yourself against Zika by wearing long pants and long
sleeve shirts? */
/*
0	not recorded	-99
114	always			4
220	almost always	5
356	sometimes		6
162	rarely			7
72	never			8
*/

recode Q392 (4=5 "Always") (5=4 "Almost always") (6=3 "Sometimes") (7=2 "Rarely") (8=1 "Never") (-99=.), gen(protect_repel_w2)
/* How often do you protect yourself against Zika by using mosquito repellent?*/
/*
3	not recorded	-99
157	always			4
166	almost always	5
263	sometimes		6
217	rarely			7
118	never			8
*/

recode Q394 (4=5 "Always") (5=4 "Almost always") (6=3 "Sometimes") (7=2 "Rarely") (8=1 "Never") (-99=.), gen(protect_screen_w2)
/* How often do you protect yourself against Zika by using screens or closing
windows? */
/*
2	not recorded	-99
227	always			4
207	almost always	5
175	sometimes		6
157	rarely			7
156	never			8
*/

recode Q398 (1=1 "Yes") (2=0 "No") (-99=.), gen(has_zika_w2)
/* Has a doctor or health professional ever told you that you have contracted
 the Zika virus? */
/*
1		not recorded	-99
113		yes				1
808		no				2
*/

recode Q402 (1=1 "Yes") (2=0 "No") (-99=.), gen(family_zika_w2)
/* To the best of your knowledge, has a doctor or health professional ever told
any of your closest friends or family members that he or she contracted the Zika
virus? */
/*
12		not recorded	-99
264		yes				1
646		no				2
*/

recode Q72 (1=1 "Yes") (2=0 "No") (-99=.), gen(outsidehelp2_w2)
/* Did you make an effort to look for information during the study? */
/*
1		not recorded	-99
209		yes				1
713		no				2
*/

recode Q74 (1=1 "Yes") (2=0 "No") (-99=.), gen(more_info_w2)
/* Indicate if you are interested in reading more information about Zika and 
how to prevent it. */
/*
1	not recorded	-99
676	yes				1
242	no				2
*/

/* Clean Wave 2 */ 

* drop if they NEVER reach the main survey section *
gen reachexp = (Q338==1)
drop if reachexp == 0

* characters before char 15 are meaningless *
gen id = substr(psid, 15, .)

* tag duplicates *
duplicates tag id, gen(duplicate)
list psid V8 if duplicate>0
codebook id
tab duplicate

* keep only the first duplicate *
bysort id: egen mintime = min(V8) /*min start time*/
format mintime %tc

/* due to Stata quirk, mintime is within 90 seconds of true minimum. Keep 
observations within the 90-second threshold*/
gen keeper=(abs((V8-mintime)/1000)<=90)
drop if keeper==0 & duplicate>0
sort id V8
list id V8 mintime keeper Q3 Q340 Q362 Q394 if duplicate>0
drop duplicate
duplicates tag id, gen(duplicate)
list id V8 mintime keeper Q3 Q340 Q362 Q394 if duplicate>0

/* drop problematic respondents who took Wave 2 survey multiple times in the 
span of a minute */ 
drop if id == "RKNvJArg**"
drop if id == "tedHw7JQ**"

/* drop problematic respondent who will be dropped from Wave 1 data */ 
drop if id == "o_vq6I4Q**"

list id V8 mintime Q33 if duplicate>0

gen w2start=V8

save "Zika_CT_experimentwave_2-merge.dta", replace

*** WAVE 1 ***

/* check wave 1 PSIDs with master list */

clear
import delimited "PSIDs for recontact.csv"
rename v1 psid
duplicates drop psid, force
save "PSIDs_for_recontact.dta", replace

clear
use "Zika_CT_experimentwave_1.dta"

/* Recode Wave 1 variables */

rename treat1 control
/* Neutral control article */
/*
667		1	received treatment
2990	.	did not receive treatment
*/
recode control (.=0)

rename treat2 mosquito
/* Facts about Aedes aegypti */
/*
631		1	received treatment
3026	.	did not receive treatment
*/
recode mosquito (.=0)

rename treat3 behavior
/* How to avoid mosquito breeding sites in your home  */
/*
647		1	received treatment
3010	.	did not receive treatment
*/
recode behavior (.=0)

rename treat4 correction
/* Facts countering Zika rumors */
/*
616		1	received treatment
3041	.	did not receive treatment
*/
recode correction (.=0)

gen cond = .
replace cond = 1 if correction == 1
replace cond = 2 if control == 1
replace cond = 3 if mosquito == 1
replace cond = 4 if behavior == 1
tab cond, missing

recode Q3 (1=1 "Yes") (2=0 "No"), gen(minor)

/* Are you 18 or older? */
/*
3140	Yes	1
445		No	2
*/

rename Q4 age

recode Q5 (1=0 "Male") (2=1 "Female") (3=2 "Other"), gen(sex)
/*
1099	1 	Male
1353	2 	Female
6		3	Other
*/

rename Q6 mother_education 
/*
47	15	none
265	16	primary, incomplete
307	17	primary, complete
187	18	secondary, incomplete
607	19	secondary, complete
36	20	technological, incomplete
160	21	technological, complete
172	22	university, incomplete
677	23	university, complete
*/

sort psid
merge m:1 psid using "ssi-education.dta"
tab _merge /*none unmatched*/
drop _merge

rename Q7 race 
/*
1519	1	white
27		2	indigenous
163		3	black
654		4	brown
68		5	yellow
27		6	other
*/
gen white = .
replace white = 0 if race != 1 & race != .
replace white = 1 if race == 1
gen black = .
replace black = 1 if race == 3
gen pardo = .
replace pardo = 1 if race == 4 

rename Q8 income
/*
18		-99	not recorded
27		18	no income
57		19	less than 700R$
96		20	700R$-950R$
153		21	951R$-1050R$
116		22	1051R$-1200R$
55		23	1201R$-1350R$
79		24	1351R$-1500R$
59		25	1501R$-1750R$
97		26	1751R$-1950R$
112		27	1951R$-2150R$
108		28	2151R$-2350R$
75		29	2351R$-2550R$
209		30	2551R$-3150R$
168		31	3151R$-3800R$
217		32	3801R$-4950R$
362		33	4951R$-6700R$
658		34	more than 6700R$
*/
replace income=. if income==-99
su income, detail

rename Q9 state
/*
4	Acre 				1
17 	Alagoas 			2
5 	Amap¬∑				3
19 	Amazonas			4
135 Bahia				5
68 	Cear¬∑				6
65 	Distrito Federal	7
58 	Esp√årito Santo		8
64 	Goi¬∑s				9
19 	Maranh‚Äûo			10
26 	Mato Grosso			11
18 	Mato Grosso do Sul	12
273 Minas Gerais		13	
39 	Par¬∑				14
28 	Para√åba			15
140 Paran¬∑				16
101 Pernambuco			17
13 	Piaui				18
271 Rio de Janeiro		19
34 	Rio Grande do Norte	20
153 Rio Grande do Sul	21
12 	Rond√ônia			22
2 	Roraima				23
97 	Santa Catarina		24
779 S‚Äûo Paulo			25
10 	Sergipe				26
8 	Tocantins			27
*/

rename Q10 postal_code

rename Q11 religion
/*
7		not recorded		-99
1057	Catholic			13
227		Protestant			14
14		other non-Christian	15
333		Pentecostal			16
18		Mormon				17
46		traditional			18
222		kardecism			19
9		Jew					20
26		Jehovah's witness	21
294		none				22
134		agnostic			23
71		other				24
*/
replace religion=. if religion==-99

gen catholic = .
replace catholic = 1 if religion == 13

gen employed = .
replace employed = 1 if Q12 == 8
replace employed = 0 if (Q12 == 9 | Q12 == 10 | Q12 == 11 | Q12 == 12 | Q12 == 13 | Q12 == 32)

rename Q12 employment
/*
10		not recorded		-99
1525	employed			8
198		unemployed with job	9
245		looking				10
223		student				11
114		homemaker			12
111		disabled			13
32		unemployed			32
*/
replace employment=. if employment==-99

recode Q13 (1=4 "Very interested") (2=3 "Somewhat interested") (3=2 "A little interested") (4=1 "Not interested") (-99=.), gen(political_interest)
/*
61	not recorded		-99
840	very interested		1
816	somewhat interested	2
590	a little interested	3
151	not interested 		4
*/

recode Q17 (1=5 "Strongly approve") (2=4 "Somewhat approve") (4=2 "Somewhat disapprove") (5=1 "Strongly disapprove") (-99=.), gen(govt_approval)
/*
51	not recorded					-99
83	strongly approve				1
319	somewhat approve				2
571	neither approve nor disapprove	3
539	somewhat disapprove				4
895	strongly disapprove				5
*/

rename Q16_1 political_beliefs
* 1 is very left-leaning, 10 is very right-leaning *
/*
150	1	1
100	2	2
123	3	3
158	4	4
752	5	5
480	6	6
233	7	7
206	8	8
95	9	9
161	10	10
*/

su political_beliefs, detail
tab political_beliefs

gen left=(political_beliefs<6) if political_beliefs!=.
tab left

gen right=(political_beliefs>5) if political_beliefs!=.
tab right

recode Q14 (2=0 "No") (-99=.), gen(party_member)
/*
581		Yes	1
1877	No	2
*/

rename Q15 party_affiliation
/*
146	PT (Partido dos Trabalhadores)					15
47	PMDB (Partido do Movimento Democr¬∑tico)		16
99	PSDB (Partido da Social Democracia Brasileira)	17
12	PSB (Partido Socialista Brasileiro)				18
24	Democratas										19
10	PC do B (Partido Comunista do Brasil)			20
2	PPS (Partido Popular Socialista)				21
6	PTB (Partido Trabalhista Brasileiro)			22
64	PSOL (Partido Socialismo e Liberdade)			23
18	PP (Partido Progressista)						24
2	PL (Partido Liberal)							25
34	PV (Partido Verde)								26
14	PDT (Partido Democr¬∑tico Trabalhista)			27
103	Other											28
*/

recode Q125 (9=1 "Not accurate") (10=2 "A little accurate") (11=3 "More or less accurate") (12=4 "Very accurate") (-99=.), gen(cabals)
/* Major events like wars/recession/election are controlled by small groups of
people who are working secretly against the rest of us / the population */
/*
3	not recorded			-99
439	not accurate			9
536	a little accurate		10
941	more or less accurate	11
539	very accurate			12
*/

recode Q126 (9=1 "Not accurate") (10=2 "A little accurate") (11=3 "More or less accurate") (12=4 "Very accurate") (-99=.), gen(esp)
/* Some people can read other people's minds or see events before they happen */
/*
1	not recorded			-99
835	not accurate			9
506	a little accurate		10
771	more or less accurate	11
345	very accurate			12
*/

recode Q127 (14=1 "Not accurate") (15=2 "A little accurate") (16=3 "More or less accurate") (17=4 "Very accurate") (-99=.), gen(spirits)
/* Spirits can communicate with the living and act in the physical world */
/*
2	not recorded			-99
790	not accurate			14
441	a little accurate		15
738	more or less accurate	16
487	very accurate			17
*/

recode Q20_4 (6=1 "Not accurate") (7=2 "A little accurate") (8=3 "More or less accurate") (9=4 "Very accurate") (-99=.), gen(solveproblems)
/* I can solve the problems I have */
/*
6		not recorded			-99
56		not accurate			6
218		a little accurate		7
1094	more or less accurate	8
1219	very accurate			9
*/

recode Q20_11 (6=1 "Not accurate") (7=2 "A little accurate") (8=3 "More or less accurate") (9=4 "Very accurate") (-99=.), gen(pushed)
/* Sometimes, I feel like I'm being pushed */
/*
15	not recorded			-99
584	not accurate			6
787	a little accurate		7
924	more or less accurate	8
283	very accurate			9
*/

recode Q20_12 (6=1 "Not accurate") (7=2 "A little accurate") (8=3 "More or less accurate") (9=4 "Very accurate") (-99=.), gen(littlecontrol)
/* I have little control over what happens to me */
/*
11	not recorded			-99
919	not accurate			6
875	a little accurate		7
611	more or less accurate	8
177	very accurate			9
*/

recode Q20_13 (6=1 "Not accurate") (7=2 "A little accurate") (8=3 "More or less accurate") (9=4 "Very accurate") (-99=.), gen(doanything)
/* I can do just about anything I set my mind to */
/*
11		not recorded			-99
194		not accurate			6
495		a little accurate		7
1175	more or less accurate	8
718		very accurate			9
*/

recode Q20_14 (6=1 "Not accurate") (7=2 "A little accurate") (8=3 "More or less accurate") (9=4 "Very accurate") (-99=.), gen(powerless)
/* Many times, I feel powerless to deal with life's problems */
/* 
13		not recorded			-99
378		not accurate			6
705		a little accurate		7
1000	more or less accurate	8
497		very accurate			9
*/

recode Q20_15 (6=1 "Not accurate") (7=2 "A little accurate") (8=3 "More or less accurate") (9=4 "Very accurate") (-99=.), gen(future)
/* What happens to me in the future depends mostly on me */
/*
8		not recorded			-99
69		not accurate			6
196		a little accurate		7
847		more or less accurate	8
1473	very accurate			9
*/

recode Q20_16 (6=1 "Not accurate") (7=2 "A little accurate") (8=3 "More or less accurate") (9=4 "Very accurate") (-99=.), gen(importantchange)
/* There is little I can do to change the important things in my life */
/*
7		not recorded			-99
1016	not accurate			6
784		a little accurate		7
570		more or less accurate	8
216		very accurate			9
*/

recode Q22 (1=1 "Yes") (2=0 "No") (-99=.), gen(outsidehelp1)
/* Will you answer the following questions without help from outside sources? */
/*
8	not recorded	-99
1	Yes				2532
2	No				51
*/

recode Q23_1 (1=1 "True") (2=0 "False") (-99=.), gen(earthcenter)
/* The center of the earth is very hot */
/*
10		not recorded	-99
2353	True		1
224		False		2
*/

recode Q23_2 (1=1 "True") (2=0 "False") (-99=.), gen(radioactivity)
/* All radioactivity is produced by man */
/*
11		not recorded	-99
515		True			1
2061	False			2
*/

recode Q23_3 (1=1 "True") (2=0 "False") (-99=.), gen(lasers)
/* Lasers work by concentrating sound waves */
/*
15		not recorded	-99
843		True			1
1729	False			2
*/

recode Q23_4 (1=1 "True") (2=0 "False") (-99=.), gen(electrons)
/* Electrons are smaller than atoms */
/*
9		not recorded	-99
1613	True			1
965		False			2
*/

recode Q23_5 (1=1 "True") (2=0 "False") (-99=.), gen(chromosomes)
/* The father's genes decide whether a baby is a boy or girl */
/*
11		not recorded	-99
1563	True			1
1013	False			2
*/

recode Q23_6 (1=1 "True") (2=0 "False") (-99=.), gen(antibiotics)
/* Antibiotics kill viruses as well as bacteria */
/*
10		not recorded	-99
1399	True			1
1178	False			2
*/

recode Q25 (1=1 "Correct") (2=0 "Incorrect") (-99=.), gen(eartharoundsun)
/* Does the earth revolve around the sun, or does the sun revolve around the 
earth? */
/*
2		not recorded				-99
2308	earth revolves around sun	1
276		sun revolves around earth	2
*/

recode Q24 (1=0 "Incorrect") (2=0 "Incorrect") (3=1 "Correct") (-99=.), gen(yearlength)
/* How long does it take for the earth to go around the sun? */
/*
18	not recorded	-99
770		one day		1
63		one month	2
1733	one year	3
*/

recode Q26 (1=0 "Incorrect") (2=0 "Incorrect") (3=0 "Incorrect") (5=1 "Correct") (-99=.), gen(uspresident)
/* What is the name of the current president of the United States of America? */
/*
3		not recorded	-99
43		Barack Obama	1
31		George W. Bush	2
5		Al Gore			3
2502	Donald Trump	5
*/

recode Q27 (1=1 "Correct") (2=0 "Incorrect") (3=0 "Incorrect") (4=0 "Correct") (-99=.), gen(nigeria)
/* On what continent is Nigeria? */
/*
2		not recorded	-99
2313	Africa			1
130		Asia			2
43		North America	3
96		South America	4
*/

recode Q28 (1=0 "Incorrect") (2=0 "Incorrect") (3=1 "Correct") (4=0 "Incorrect") (-99=.), gen(presidentialterm)
/* The presidential term in Brazil is how long? */
/*
2		not recorded	-99
20		Two years		1
24		Three years		2
2474	Four years		3
64		Five years		4
*/

recode Q29 (1=0 "Incorrect") (2=0 "Incorrect") (3=1 "Correct") (4=0 "Incorrect") (-99=.), gen(howmanyreps)
/* How many representatives are in the Chamber of Deputies? */
/*
7		not recorded	-99
277		81				1
473		232				2
371		388				3
1456	513				4
*/

recode Q30_1 (1=1 "Not at all") (2=2 "Somewhat closely") (3=3 "Closely") (4=4 "Very closely"), gen(news_politics)
/* How closely do you follow political news? */
/*
14		not recorded		-99
173		Not at all			1	
1067	Somewhat closely	2
864		Closely				3
462		Very closely		4
*/
recode news_politics (-99=.)

recode Q30_2 (1=1 "Not at all") (2=2 "Somewhat closely") (3=3 "Closely") (4=4 "Very closely"), gen(news_sports)
/* How closely do you follow sports news? */
/*
15		not recorded		-99
620		Not at all			1	
1001	Somewhat closely	2
504		Closely				3
440		Very closely		4
*/
recode news_sports (-99=.)

recode Q30_3 (1=1 "Not at all") (2=2 "Somewhat closely") (3=3 "Closely") (4=4 "Very closely"), gen(news_religion)
/* How closely do you follow religious news? */
/*
12		not recorded		-99
515		Not at all			1	
1213	Somewhat closely	2
545		Closely				3
295		Very closely		4
*/
recode news_religion (-99=.)

recode Q30_4 (1=1 "Not at all") (2=2 "Somewhat closely") (3=3 "Closely") (4=4 "Very closely"), gen(news_world)
/* How closely do you follow foreign affairs news? */
/*
15		not recorded		-99
136		Not at all			1	
1232	Somewhat closely	2
874		Closely				3
323		Very closely		4
*/
recode news_world (-99=.)

recode Q30_5 (1=1 "Not at all") (2=2 "Somewhat closely") (3=3 "Closely") (4=4 "Very closely"), gen(news_science)
/* How closely do you follow science news? */
/*
11		not recorded		-99
230		Not at all			1	
1307	Somewhat closely	2
739		Closely				3
293		Very closely		4
*/
recode news_science (-99=.)

recode Q30_6 (1=1 "Not at all") (2=2 "Somewhat closely") (3=3 "Closely") (4=4 "Very closely"), gen(news_tech)
/* How closely do you follow technology news? */
/*
14		not recorded		-99
119		Not at all			1	
1103	Somewhat closely	2
887		Closely				3
457		Very closely		4
*/
recode news_tech (-99=.)

recode Q30_7 (1=1 "Not at all") (2=2 "Somewhat closely") (3=3 "Closely") (4=4 "Very closely"), gen(news_entertainment)
/* How closely do you follow entertainment news? */
/*
12		not recorded		-99
656		Not at all			1	
1275	Somewhat closely	2
443		Closely				3
194		Very closely		4
*/
recode news_entertainment (-99=.)

recode Q31_1 (5=1 "Not interested") (6=2 "Slightly interested") (7=3 "More than slightly interested") (8=4 "Very interested") (-99=.), gen(interest_politics)
/* How interested are you in government and politics? */
/*
62	not recorded		-99
281	Not interested		5
778	Slightly interested	6
853	Somewhat interested	7
602	Very interested		8
*/

recode Q31_2 (5=1 "Not interested") (6=2 "Slightly interested") (7=3 "More than slightly interested") (8=4 "Very interested") (-99=.), gen(interest_sports)
/* How interested are you in sports? */
/*
59	not recorded		-99
612	Not interested		5
831	Slightly interested	6
541	Somewhat interested	7
533	Very interested		8
*/

recode Q31_3 (5=1 "Not interested") (6=2 "Slightly interested") (7=3 "More than slightly interested") (8=4 "Very interested") (-99=.), gen(interest_religion)
/* How interested are you in religion? */
/*
63	not recorded		-99
449	Not interested		5
945	Slightly interested	6
641	Somewhat interested	7
478	Very interested		8
*/

recode Q31_4 (5=1 "Not interested") (6=2 "Slightly interested") (7=3 "More than slightly interested") (8=4 "Very interested") (-99=.), gen(interest_world)
/* How interested are you in international affairs? */
/*
68	not recorded		-99
134	Not interested		5
885	Slightly interested	6
984	Somewhat interested	7
505	Very interested		8
*/

recode Q31_5 (5=1 "Not interested") (6=2 "Slightly interested") (7=3 "More than slightly interested") (8=4 "Very interested") (-99=.), gen(interest_science)
/* How interested are you in scientific discoveries? */
/* 
63	not recorded		-99
117	Not interested		5
756	Slightly interested	6
952	Somewhat interested	7
688	Very interested		8
*/

recode Q31_6 (5=1 "Not interested") (6=2 "Slightly interested") (7=3 "More than slightly interested") (8=4 "Very interested") (-99=.), gen(interest_tech)
/* How interested are you in technology? */
/*
59	not recorded		-99
71	Not interested		5
679	Slightly interested	6
909	Somewhat interested	7
858	Very interested		8
*/

recode Q31_7 (5=1 "Not interested") (6=2 "Slightly interested") (7=3 "More than slightly interested") (8=4 "Very interested") (-99=.), gen(interest_entertainment)
/* How interested are you in entertainment or celebrities? */
/*
61		not recorded		-99
766		Not interested		5
1065	Slightly interested	6
450		Somewhat interested	7
234		Very interested		8
*/

rename Q136 confidence_court 
/* How much confidence do you have in the Supreme Federal Court? */
/*
19	not recorded	-99
470	No confidence	1
428					2
428					3
489					4
418					5
208					6
110	Full confidence	7
*/

rename Q137 confidence_congress
/* How much confidence do you have in the National Congress? */
/*
16		not recorded	-99
1221	No confidence	1
540						2
347						3
255						4
113						5
43						6
35		Full confidence	7
*/

rename Q138 confidence_press
/* How much confidence do you have in the press? */
/*
11		not recorded	-99
365		No confidence	1
434						2
524						3
578						4
404						5
156						6
98		Full confidence	7
*/

rename Q139 confidence_govt
/* How much confidence do you have in the government? */
/*
20		not recorded	-99
1199	No confidence	1
580						2
338						3
264						4
107						5
31						6
31		Full confidence	7
*/

rename Q140 confidence_medicine
/* How much confidence do you have in medicines [doctors? portuguese word is "remedios"] ? */
/*
17		not recorded	-99
44		No confidence	1
116						2
277						3
549						4
749						5
554						6
264		Full confidence	7
*/

rename Q141 confidence_scientists
/* How much confidence do you have in scientists? */
/*
19		not recorded	-99
40		No confidence	1
102						2
168						3
396						4
718						5
705						6
422		Full confidence	7
*/

rename Q142 confidence_banks
/* How much confidence do you have in banks and financial institutions? */
/*
20		not recorded	-99
309		No confidence	1
400						2
489						3
604						4
459						5
209						6
80		Full confidence	7
*/

rename Q143 confidence_govhealth
/* How much confidence do you have in the people who manage the Ministry of 
Health */
/*
12		not recorded	-99
602		No confidence	1
564						2
514						3
474						4
260						5
95						6
49		Full confidence	7
*/

foreach var of varlist confidence* {
replace `var'=. if `var'==-99
}

recode Q43 (6=1 "Not accurate") (10=2 "Somewhat accurate") (11=3 "More or less accurate") (12=4 "Very accurate") (-99=.), gen(zika_bite)
/* A person may be infected with Zika by the bite of a mosquito carrying the
virus. */
/*
13		not recorded			-99
48		Not accurate			6
64		Somewhat accurate		10
176		More or less accurate	11
2247	Very accurate			12
*/

recode Q44 (10=1 "Not accurate") (11=2 "Somewhat accurate") (12=3 "More or less accurate") (13=4 "Very accurate") (-99=.), gen(zika_neuro)
/* Being infected with the Zika virus increases the likelihood of developing 
neurological problems. */
/*
13		not recorded			-99
322		Not accurate			10
430		Somewhat accurate		11
786		More or less accurate	12
997		Very accurate			13
*/

recode Q45 (10=1 "Not accurate") (11=2 "Somewhat accurate") (12=3 "More or less accurate") (13=4 "Very accurate") (-99=.), gen(zika_sex)
/* A person can become infected with Zika by having sex with someone who is
infected. */
/*
15		not recorded			-99
1405	Not accurate			10
316		Somewhat accurate		11
405		More or less accurate	12
407		Very accurate			13
*/

recode Q46 (10=4 "Not accurate") (11=3 "Somewhat accurate") (12=2 "More or less accurate") (13=1 "Very accurate") (-99=.), gen(zika_handshake)
/* A person can become infected with Zika through casual contact with someone 
who is infected.
*/
/*
10		not recorded			-99
2227	Not accurate			10
160		Somewhat accurate		11
92		More or less accurate	12
59		Very accurate			13
*/

recode Q47 (10=1 "Not accurate") (11=2 "Somewhat accurate") (12=3 "More or less accurate") (13=4 "Very accurate") (-99=.), gen(zika_gmo)
/* Genetically modified mosquitoes - not wild mosquitoes - caused Zika's
 outbreak in Latin America. */
/* 
14		not recorded			-99
1227	Not accurate			10
511		Somewhat accurate		11
498		More or less accurate	12
298		Very accurate			13
*/

recode Q49 (10=1 "Not accurate") (11=2 "Somewhat accurate") (12=3 "More or less accurate") (13=4 "Very accurate") (-99=.), gen(micro_larvicide)
/* The prevalence of microcephaly in Brazil has increased as a result of 
larvicides added to water to kill mosquitoes */
/*
9		not recorded			-99
1574	Not accurate			10
430		Somewhat accurate		11
337		More or less accurate	12
173		Very accurate			13
*/

recode Q50 (10=1 "Not accurate") (11=2 "Somewhat accurate") (12=3 "More or less accurate") (13=4 "Very accurate") (-99=.), gen(micro_vaccine)
/* The prevalence of microcephaly in Brazil has increased as a result of the
government determining increased use of vaccines among pregnant women. */
/*
10		not recorded			-99
1726	Not accurate			10
348		Somewhat accurate		11
252		More or less accurate	12
187		Very accurate			13
*/

recode Q53 (4=1) (5=2) (6=3) (7=4) (8=5) (9=6) (10=7) (11=8) (12=9) (13=10) (-99=.), gen(policy_dtap)
/* Do you approve or disapprove of the government's recommendation that pregnant
women receive the DTaP vaccine? */
/* 1 = strongly disapprove, 10 = strongly approve */
/*
20	not recorded		-99
77	strongly disapprove	1
43						2
64						3
95						4
358						5
228						6
263						7
378						8
246						9
730	strongly approve	10
*/

recode Q133 (4=1) (5=2) (6=3) (7=4) (8=5) (9=6) (10=7) (11=8) (12=9) (13=10) (-99=.), gen(policy_larvicide)
/* Do you approve or disapprove of the use of larvicides in water to prevent 
mosquito breeding in areas where Zika is believed to be present? */
/* 1 = strongly disapprove, 10 = strongly approve */
/*
13	not recorded		-99
97	strongly disapprove	1
48						2
81						3
117						4
273						5
213						6
279						7
323						8
280						9
793	strongly approve	10
*/

recode Q134 (4=1) (5=2) (6=3) (7=4) (8=5) (9=6) (10=7) (11=8) (12=9) (13=10) (-99=.), gen(policy_gmo)
/* Do you approve or disapprove of the release of genetically modified 
mosquitoes to the environment to combat the spread of Zika? */
/* 1 = strongly disapprove, 10 = strongly approve */
/*
13	not recorded		-99
173	strongly disapprove	1
77						2
86						3
91						4
289						5
229						6
267						7
335						8
239						9
718	strongly approve	10
*/

recode Q135 (4=1) (5=2) (6=3) (7=4) (8=5) (9=6) (10=7) (11=8) (12=9) (13=10) (-99=.), gen(policy_homesearch)
/* Do you approve or disapprove of government health workers entering private 
homes to look for mosquito breeding points / outbreaks when the resident or 
owner is not present? */
/* 1 = strongly disapprove, 10 = strongly approve */
/*
12	not recorded		-99
234	strongly disapprove	1
71						2
85						3
78						4
209						5
176						6
199						7
300						8
245						9
908	strongly approve	10
*/

recode Q64 (7=1 "No risk") (8=2 "Very low risk") (9=3 "Low risk") (10=4 "Moderate risk") (12=5 "High risk") (13=6 "Very high risk") (-99=.), gen(risk_human)
/* How much risk do you think the Zika virus presents to human health in
Brazil? */
/*
5	not recorded	-99
16	no risk			7
33	very low risk	8
57	low risk		9
434	moderate risk	10
982	high risk		12
990	very high risk	13
*/

recode Q65 (7=1 "No risk") (8=2 "Very low risk") (9=3 "Low risk") (10=4 "Moderate risk") (12=5 "High risk") (13=6 "Very high risk") (-99=.), gen(risk_personal)
/* How much risk do you believe the Zika virus represents to you personally? */
/*
6	not recorded	-99
68	no risk			7
227	very low risk	8
377	low risk		9
850	moderate risk	10
546	high risk		12
440	very high risk	13
*/

recode Q66 (1=1 "Not at all concerned") (2=2 "Not very concerned") (3=3 "Somewhat concerned") (4=4 "Very concerned") (5=5 "Extremely concerned") (-99=.), gen(concern)
/* How concerned are you personally about the potential health risk for Zika? */
/* 
36	not recorded			-99
88	not at all concerned	1
281	not very concerned		2
774	somewhat concerned		3
821	very concerned			4
511	extremely concerned		5
*/

recode Q128 (4=5 "Always") (5=4 "Almost always") (6=3 "Sometimes") (7=2 "Rarely") (8=1 "Never") (-99=.), gen(protect_sleeves)
/* How often do you protect yourself against Zika by wearing long pants and long
sleeve shirts? */
/*
8	not recorded	-99
279	always			4
549	almost always	5
899	sometimes		6
521	rarely			7
253	never			8
*/

recode Q129 (4=5 "Always") (5=4 "Almost always") (6=3 "Sometimes") (7=2 "Rarely") (8=1 "Never") (-99=.), gen(protect_repel)
/* How often do you protect yourself against Zika by using mosquito repellent?*/
/*
22	not recorded	-99
395	always			4
466	almost always	5
671	sometimes		6
574	rarely			7
381	never			8
*/

recode Q130 (4=5 "Always") (5=4 "Almost always") (6=3 "Sometimes") (7=2 "Rarely") (8=1 "Never") (-99=.), gen(protect_screen)
/* How often do you protect yourself against Zika by using screens or closing
windows? */
/*
12	not recorded	-99
622	always			4
542	almost always	5
482	sometimes		6
406	rarely			7
445	never			8
*/

recode Q69 (1=1 "Yes") (2=0 "No") (-99=.), gen(has_zika)
/* Has a doctor or health professional ever told you that you have contracted
 the Zika virus? */
/*
6		not recorded	-99
290		yes				1
2211	no				2
*/

recode Q70 (1=1 "Yes") (2=0 "No") (-99=.), gen(family_zika)
/* To the best of your knowledge, has a doctor or health professional ever told
any of your closest friends or family members that he or she contracted the Zika
virus? */
/*
48		not recorded	-99
665		yes				1
1789	no				2
*/

recode Q72 (1=1 "Yes") (2=0 "No") (-99=.), gen(outsidehelp2)
/* Did you make an effort to look for information during the study? */
/*
1		not recorded	-99
207		yes				1
1035	no				2
*/

recode Q74 (1=1 "Yes") (2=0 "No") (-99=.), gen(moreinfo)
/* Indicate if you are interested in reading more information about Zika and 
how to prevent it. */
/*
9	not recorded	-99
787	yes				1
443	no				2
*/

/* Clean Wave 1 */

count

* drop if they NEVER reach the experimental manipulation *
gen reachexp = (Q33==1)
drop if reachexp == 0

* characters before char 15 are meaningless *
gen id = substr(psid, 15, .)

/* drop duplicates of id, keeping first finished observation in each group */

count

* tag duplicates *
duplicates tag id, gen(duplicate)
list id V8 if duplicate>0
codebook id
tab duplicate
duplicates list id

* keep only the first duplicate *
bysort id: egen mintime = min(V8)
format mintime %tc
gen keeper=(abs((V8-mintime)/1000)<=90)
list id V8 mintime keeper if duplicate>0
drop if keeper==0 & duplicate>0
list id V8 mintime Q33 if duplicate>0
codebook id
tab duplicate
duplicates list id

/* drop problematic respondent who took survey multiple times with sessions
beginning just seconds apart */ 
drop if id == "o_vq6I4Q**"
drop if id == "RKNvJArg**"
drop if id == "tedHw7JQ**"

count

drop if cond==.

count

* merge wave 1 and master datasets
merge 1:1 psid using PSIDs_for_recontact.dta

/* drop same problematic respondent in merged data */
drop if _merge == 2

gen invited = (_merge == 3)
drop _merge

save "Zika_CT_experimentwave_1-merge.dta", replace

/* merge wave 1 with wave 2 */
clear
use "Zika_CT_experimentwave_1-merge.dta"
duplicates list id

merge 1:1 id using "Zika_CT_experimentwave_2-merge"
tab _merge
list id V8 if _merge==2

/*remove wave 2 responses from mistakenly invited participants*/
foreach var of varlist *_w2 {
replace `var'=. if psid=="00UPWNOQI_HBawDM4dxtPA**" | psid=="00UPWNOQI_EbXbOAksNcLw**" 
}

tab correction control if zika_sex_w2!=. | zika_gmo_w2!=., missing

duplicates report psid

save "Zika_CT_experiment-merged.dta", replace

*Table S2
tab Q5
tab white
tab age
tab education, missing
gen hs_degree=(strpos(education,"Highschool complete")!=0 | strpos(education,"University complete")!=0)
tab hs_degree if educ!=""
tab income
tab political_beliefs
su cabal

*response time
gen response_time=.
replace response_time=Q41_3 if correction==1
replace response_time=Q35_3 if correction==0 & mosquito==0 & behavior==0
xtile response_quartile_exp=response_time,nq(4)

egen prior_response_time=rowtotal(Q108_3 Q106_3 Q109_3 Q110_3 Q111_3 Q118_3 Q119_3 Q112_3 Q147_3 Q114_3 Q115_3 Q116_3 Q117_3)
su prior_response_time, detail

gen very_fast=prior<193.629 /*10th percentile fastest*/
xtile response_quartile=prior_response_time,nq(4)

/** Creating scales needed in analysis **/
/*note: variables collected in wave 2 end in "_w2"; otherwise from wave 1*/

* Zika knowledge (waves 1 and 2):
factor zika_bite zika_neuro zika_sex zika_handshake, pcf 
rotate, varimax 
alpha zika_bite zika_neuro zika_sex zika_handshake, item casewise gen(ZikaKnow) asis 

factor zika_bite_w2 zika_neuro_w2 zika_sex_w2 zika_handshake_w2, pcf  
rotate, varimax
alpha zika_bite_w2 zika_neuro_w2 zika_sex_w2 zika_handshake_w2, item casewise gen(ZikaKnow_w2) asis 

* Zika misperceptions:
factor zika_gmo micro_larvicide micro_vaccine, pcf
rotate, varimax 
alpha zika_gmo micro_larvicide micro_vaccine, item casewise gen(ZikaMisp) 

factor zika_gmo_w2 micro_larvicide_w2 micro_vaccine_w2, pcf
rotate, varimax 
alpha zika_gmo_w2 micro_larvicide_w2 micro_vaccine_w2, item casewise gen(ZikaMisp_w2) 

*Policy opinions: 
factor policy_gmo policy_dtap policy_larvicide, pcf 
rotate, varimax 
alpha policy_gmo policy_dtap policy_larvicide, item casewise gen(ZikaOpinions) 

factor policy_gmo_w2 policy_dtap_w2 policy_larvicide_w2, pcf 
rotate, varimax 
alpha policy_gmo_w2 policy_dtap_w2 policy_larvicide_w2, item casewise gen(ZikaOpinions_w2) 

*Behavioral intentions: 
factor protect_sleeves protect_repel protect_screen, pcf 
rotate, varimax 
alpha protect_sleeves protect_repel protect_screen, item casewise gen(ZikaBehavior) 

factor protect_sleeves_w2 protect_repel_w2 protect_screen_w2, pcf 
rotate, varimax 
alpha protect_sleeves_w2 protect_repel_w2 protect_screen_w2, item casewise gen(ZikaBehavior_w2) 

*Science and political knowledge (latter not used)
gen sciknow1=(Q23_1==1)
gen sciknow2=(Q23_2==2)
gen sciknow3=(Q23_3==2)
gen sciknow4=(Q23_4==1)
gen sciknow5=(Q23_5==1)
gen sciknow6=(Q23_6==2)
gen sciknow7=(Q25==1)
gen sciknow8=(Q24==3)

gen polknow1=(Q26==5)
gen polknow2=(Q27==1)
gen polknow3=(Q28==3)
gen polknow4=(Q29==3)

factor sciknow* polknow*, pcf
rotate, varimax

alpha sciknow*, item 
alpha polknow*, item  
egen sciknow=rowtotal(sciknow*)
replace sciknow=sciknow/8
egen polknow=rowtotal(polknow*)
replace polknow=polknow/4

/** ANALYSIS **/

reg ZikaMisp correction if correction==1 | control==1, robust 
est store S1M

reg zika_gmo correction if correction==1 | control==1, robust 
est store S1A

reg micro_larvicide correction if correction==1 | control==1, robust 
est store S1B

reg micro_vaccine correction if correction==1 | control==1, robust 
est store S1C

*Table 2a
estout S1M S1A S1B S1C, style(tex) replace varwidth(25) collabels("") cells(b(star fmt(%9.2f)) se(par fmt(%9.2f))) stats(N, fmt(%9.2f %9.0f) labels("N")) starlevels(* 0.05 ** 0.01 *** 0.005) 

oprobit zika_gmo correction if correction==1 | control==1, robust 
est store S1A

oprobit micro_larvicide correction if correction==1 | control==1, robust 
est store S1B

oprobit micro_vaccine correction if correction==1 | control==1, robust 
est store S1C

*Table S5a
estout S1A S1B S1C, style(tex) replace varwidth(25) collabels("") cells(b(star fmt(%9.2f)) se(par fmt(%9.2f))) stats(N, fmt(%9.2f %9.0f) labels("N")) starlevels(* 0.05 ** 0.01 *** 0.005) 

reg ZikaMisp correction mosquito behavior, robust
est store S1M
lincom correction-mosquito
lincom correction-behavior
lincom mosquito-behavior

reg zika_gmo correction mosquito behavior, robust 
est store S1A
lincom correction-mosquito
lincom correction-behavior
lincom mosquito-behavior

reg micro_larvicide correction mosquito behavior, robust
est store S1B
lincom correction-mosquito
lincom correction-behavior
lincom mosquito-behavior

reg micro_vaccine correction mosquito behavior, robust 
est store S1C
lincom correction-mosquito
lincom correction-behavior
lincom mosquito-behavior

*Table S3a
estout S1M S1A S1B S1C, style(tex) replace varwidth(25) collabels("") cells(b(star fmt(%9.2f)) se(par fmt(%9.2f))) stats(N, fmt(%9.2f %9.0f) labels("N")) starlevels(* 0.05 ** 0.01 *** 0.005) 

*response times
bysort correction: su response_time if mosquito==0 & behavior==0, detail

gen cond2 = .
replace cond2 = 1 if control == 1
replace cond2 = 2 if correction == 1
replace cond2 = 3 if mosquito == 1
replace cond2 = 4 if behavior == 1

label def cond2label 1 "Control condition" 2 "Myth correction" 3 "Mosquito information" 4 "Preventive behavior"
label val cond2 cond2label

preserve
gen belief1=ZikaMisp
gen belief2=zika_gmo
gen belief3=micro_larvicide
gen belief4=micro_vaccine

reshape long belief,i(psid) j(dv)

drop if dv==1

label def dvlabel 2 "GMO mosquitoes" 3 "Larvicide harm" 4 "Vaccine harm"
label val dv dvlabel

*Figure S1a
cibar belief,  over1(cond2) over2(dv)  bargap(25) barcol(gs12 teal sand purple) gap(150) ciopts(lcolor(black)) graphopts( ytitle("") scheme(plotplainblind) ylabel(1 "Not at all accurate" 2 "Not very accurate" 3 "Somewhat accurate",angle(0) labsize(*1)) yscale(r(1 3.1)))

restore

gen zika_handshake_belief=zika_handshake
recode zika_handshake_belief (1=4) (2=3) (3=2) (4=1)

reg zika_neuro correction if correction==1 | control==1, robust 
est store B1

reg zika_bite correction if correction==1 | control==1, robust 
est store C1

reg zika_sex correction if correction==1 | control==1, robust 
est store D1

reg zika_handshake_belief correction if correction==1 | control==1, robust 
est store E1

*Table 3a
estout B1 C1 D1 E1, style(tex) replace varwidth(25) collabels("") cells(b(star fmt(%9.2f)) se(par fmt(%9.2f))) stats(N, fmt(%9.2f %9.0f) labels("N")) starlevels(* 0.05 ** 0.01 *** 0.005) 

oprobit zika_neuro correction if correction==1 | control==1, robust 
est store B1

oprobit zika_bite correction if correction==1 | control==1, robust 
est store C1

oprobit zika_sex correction if correction==1 | control==1, robust 
est store D1

oprobit zika_handshake_belief correction if correction==1 | control==1, robust 
est store E1

*Table S5b
estout B1 C1 D1 E1, style(tex) replace varwidth(25) collabels("") cells(b(star fmt(%9.2f)) se(par fmt(%9.2f))) stats(N, fmt(%9.2f %9.0f) labels("N")) starlevels(* 0.05 ** 0.01 *** 0.005) 

reg zika_neuro correction mosquito behavior, robust 
est store B1
lincom correction-mosquito
lincom correction-behavior
lincom mosquito-behavior

reg zika_bite correction mosquito behavior, robust
est store C1
lincom correction-mosquito
lincom correction-behavior
lincom mosquito-behavior

reg zika_sex correction mosquito behavior, robust 
est store D1
lincom correction-mosquito
lincom correction-behavior
lincom mosquito-behavior

reg zika_handshake_belief correction mosquito behavior, robust 
est store E1
lincom correction-mosquito
lincom correction-behavior
return list
lincom mosquito-behavior

*Table S3b
estout B1 C1 D1 E1, style(tex) replace varwidth(25) collabels("") cells(b(star fmt(%9.2f)) se(par fmt(%9.2f))) stats(N, fmt(%9.2f %9.0f) labels("N")) starlevels(* 0.05 ** 0.01 *** 0.005) 

*effect sizes
foreach var of varlist zika_neuro zika_bite zika_sex zika_handshake_belief {
esize twosample `var' if control==1 | correction==1, by(correction) cohensd 
}

gen handshake2=zika_handshake
recode handshake2 (1=4)(2=3)(3=2)(4=1)

egen averageacc=rowmean(zika_neuro zika_bite zika_sex handshake2)
reg averageacc correction mosquito behavior, robust 
esize twosample averageacc if control==1 | correction==1, by(correction) cohensd 

*joint test
mvreg zika_neuro zika_bite zika_sex zika_handshake_belief = correction if correction==1 | control==1
test correction

preserve
gen belief1=zika_neuro
gen belief2=zika_bite
gen belief3=zika_sex
gen belief4=zika_handshake

reshape long belief,i(psid) j(dv)

label def dvlabel 1 "Neurological" 2 "Mosquito" 3 "Sexual" 4 "Not casual"
label val dv dvlabel

*Figure S1b
cibar belief,  over1(cond2) over2(dv) blsize(large) bargap(25) barcol(gs12 teal sand purple) gap(150) ciopts(lcolor(black)) graphopts( ytitle("") scheme(plotplainblind) ylabel(1 "Not at all accurate" 2 "Not very accurate" 3 "Somewhat accurate" 4 "Very accurate",angle(0) labsize(*1)) yscale(r(1 4.1)))

restore

reg ZikaOpinions correction mosquito behavior, robust
est store F1
lincom correction-mosquito
lincom correction-behavior
lincom mosquito-behavior

reg ZikaBehavior correction mosquito behavior, robust
est store G1
lincom correction-mosquito
lincom correction-behavior
lincom mosquito-behavior

*Table S4
estout F1 G1, style(tex) replace varwidth(25) collabels("") cells(b(star fmt(%9.2f)) se(par fmt(%9.2f))) stats(N, fmt(%9.2f %9.0f) labels("N")) starlevels(* 0.05 ** 0.01 *** 0.005) 

gen attrition=(_merge==1) if control==1 | correction==1

*attrition analysis
ttest attrition, by(control) unequal

ttest age, by(attrition) unequal

recode sex (1=1) (2=0) (3=0) (-99=.), gen(male)
ttest male, by(attrition) unequal

ttest white, by(attrition) unequal

tab education attrition, chi2

gen w1startdate=dofc(V8)
gen w2startdate=dofc(w2start)
gen interval=w2startdate-w1startdate

*HTEs

gen correctionXcabal = correction*cabal

reg zika_neuro correction cabal correctionXcabal mosquito behavior, robust
est store B1

reg zika_bite correction cabal correctionXcabal mosquito behavior, robust
est store C1

reg zika_sex correction cabal correctionXcabal mosquito behavior, robust
est store D1

reg zika_handshake_belief correction cabal correctionXcabal mosquito behavior, robust
est store E1

reg ZikaMisp correction cabal correctionXcabal mosquito behavior, robust
est store F1

reg ZikaOpinions correction cabal correctionXcabal mosquito behavior, robust
est store G1

reg ZikaBehavior correction cabal correctionXcabal mosquito behavior, robust
est store H1

*Table S6
estout B1 C1 D1 E1 F1 G1 H1, style(tex) replace varwidth(25) collabels("") cells(b(star fmt(%9.2f)) se(par fmt(%9.2f))) stats(N, fmt(%9.0f) labels("N")) starlevels(* 0.05 ** 0.01 *** 0.005) 
 
gen correctionXconfidence_govt=correction*confidence_govt

reg zika_neuro correction confidence_govt correctionXconfidence_govt mosquito behavior, robust
est store B1

reg zika_bite correction confidence_govt correctionXconfidence_govt mosquito behavior, robust
est store C1

reg zika_sex correction confidence_govt correctionXconfidence_govt mosquito behavior, robust
est store D1

reg zika_handshake_belief correction confidence_govt correctionXconfidence_govt mosquito behavior, robust
est store E1

reg ZikaMisp correction confidence_govt correctionXconfidence_govt mosquito behavior, robust
est store F1

reg ZikaOpinions correction confidence_govt correctionXconfidence_govt mosquito behavior, robust
est store G1

reg ZikaBehavior correction confidence_govt correctionXconfidence_govt mosquito behavior, robust
est store H1

*Table S7
estout B1 C1 D1 E1 F1 G1 H1, style(tex) replace varwidth(25) collabels("") cells(b(star fmt(%9.2f)) se(par fmt(%9.2f))) stats(N, fmt(%9.0f) labels("N")) starlevels(* 0.05 ** 0.01 *** 0.005) 

gen correctionXconfidence_govhealth=correction*confidence_govhealth

reg zika_neuro correction confidence_govhealth correctionXconfidence_govhealth mosquito behavior, robust
est store B1

reg zika_bite correction confidence_govhealth correctionXconfidence_govhealth mosquito behavior, robust
est store C1

reg zika_sex correction confidence_govhealth correctionXconfidence_govhealth mosquito behavior, robust
est store D1

reg zika_handshake_belief correction confidence_govhealth correctionXconfidence_govhealth mosquito behavior, robust
est store E1

reg ZikaMisp correction confidence_govhealth correctionXconfidence_govhealth mosquito behavior, robust
est store F1

reg ZikaOpinions correction confidence_govhealth correctionXconfidence_govhealth mosquito behavior, robust
est store G1

reg ZikaBehavior correction confidence_govhealth correctionXconfidence_govhealth mosquito behavior, robust
est store H1

*Table S8
estout B1 C1 D1 E1 F1 G1 H1, style(tex) replace varwidth(25) collabels("") cells(b(star fmt(%9.2f)) se(par fmt(%9.2f))) stats(N, fmt(%9.0f) labels("N")) starlevels(* 0.05 ** 0.01 *** 0.005) 

gen correctionXconfidence_medicine=correction*confidence_medicine

reg zika_neuro correction confidence_medicine correctionXconfidence_medicine mosquito behavior, robust
est store B1

reg zika_bite correction confidence_medicine correctionXconfidence_medicine mosquito behavior, robust
est store C1

reg zika_sex correction confidence_medicine correctionXconfidence_medicine mosquito behavior, robust
est store D1

reg zika_handshake_belief correction confidence_medicine correctionXconfidence_medicine mosquito behavior, robust
est store E1

reg ZikaMisp correction confidence_medicine correctionXconfidence_medicine mosquito behavior, robust
est store F1

reg ZikaOpinions correction confidence_medicine correctionXconfidence_medicine mosquito behavior, robust
est store G1

reg ZikaBehavior correction confidence_medicine correctionXconfidence_medicine mosquito behavior, robust
est store H1

*Table S9
estout B1 C1 D1 E1 F1 G1 H1, style(tex) replace varwidth(25) collabels("") cells(b(star fmt(%9.2f)) se(par fmt(%9.2f))) stats(N, fmt(%9.0f) labels("N")) starlevels(* 0.05 ** 0.01 *** 0.005) 

gen correctionXconfidence_scientists=correction*confidence_scientists

reg zika_neuro correction confidence_scientists correctionXconfidence_scientists mosquito behavior, robust
est store B1

reg zika_bite correction confidence_scientists correctionXconfidence_scientists mosquito behavior, robust
est store C1

reg zika_sex correction confidence_scientists correctionXconfidence_scientists mosquito behavior, robust
est store D1

reg zika_handshake_belief correction confidence_scientists correctionXconfidence_scientists mosquito behavior, robust
est store E1

reg ZikaMisp correction confidence_scientists correctionXconfidence_scientists mosquito behavior, robust
est store F1

reg ZikaOpinions correction confidence_scientists correctionXconfidence_scientists mosquito behavior, robust
est store G1

reg ZikaBehavior correction confidence_scientists correctionXconfidence_scientists mosquito behavior, robust
est store H1

*Table S10
estout B1 C1 D1 E1 F1 G1 H1, style(tex) replace varwidth(25) collabels("") cells(b(star fmt(%9.2f)) se(par fmt(%9.2f))) stats(N, fmt(%9.0f) labels("N")) starlevels(* 0.05 ** 0.01 *** 0.005) 

gen college=.
replace college=0 if educationlevelbr!=""
replace college=1 if strpos(educationlevelbr,"University complete")!=0

gen correctionXcollege=correction*college

reg zika_neuro correction college correctionXcollege if correction==1 | control==1, robust
est store B1

reg zika_bite correction college correctionXcollege if correction==1 | control==1, robust
est store C1

reg zika_sex correction college correctionXcollege if correction==1 | control==1, robust
est store D1

reg zika_handshake_belief correction college correctionXcollege if correction==1 | control==1, robust
est store E1

*Table S11a
estout B1 C1 D1 E1, style(tex) replace varwidth(25) collabels("") cells(b(star fmt(%9.2f)) se(par fmt(%9.2f))) stats(N, fmt(%9.0f) labels("N")) starlevels(* 0.05 ** 0.01 *** 0.005) 

gen correctionXsciknow=correction*sciknow

reg zika_neuro correction sciknow correctionXsciknow if correction==1 | control==1, robust
est store B1

reg zika_bite correction sciknow correctionXsciknow if correction==1 | control==1, robust
est store C1

reg zika_sex correction sciknow correctionXsciknow if correction==1 | control==1, robust
est store D1

reg zika_handshake_belief correction sciknow correctionXsciknow if correction==1 | control==1, robust
est store E1

*Table S11b
estout B1 C1 D1 E1, style(tex) replace varwidth(25) collabels("") cells(b(star fmt(%9.2f)) se(par fmt(%9.2f))) stats(N, fmt(%9.0f) labels("N")) starlevels(* 0.05 ** 0.01 *** 0.005) 

reg zika_neuro_w2 correction if correction==1 | control==1, robust
est store B1

reg zika_bite_w2 correction if correction==1 | control==1, robust
est store C1

reg zika_sex_w2 correction if correction==1 | control==1, robust
est store D1

reg zika_handshake_w2 correction if correction==1 | control==1, robust
est store E1

reg ZikaMisp_w2 correction if correction==1 | control==1, robust
est store F1

reg ZikaOpinions_w2 correction if correction==1 | control==1, robust
est store G1

reg ZikaBehavior_w2 correction if correction==1 | control==1, robust
est store H1

*Table S12
estout B1 C1 D1 E1 F1 G1 H1, style(tex) replace varwidth(25) collabels("") cells(b(star fmt(%9.2f)) se(par fmt(%9.2f))) stats(N, fmt(%9.0f) labels("N")) starlevels(* 0.05 ** 0.01 *** 0.005) 

reg zika_neuro correction##ib4.response_quartile if mosquito==0 & behavior==0, robust
est store B1

reg zika_bite correction##ib4.response_quartile if mosquito==0 & behavior==0, robust
est store C1

reg zika_sex correction##ib4.response_quartile if mosquito==0 & behavior==0, robust
est store D1

reg zika_handshake_belief correction##ib4.response_quartile if mosquito==0 & behavior==0, robust
est store E1

*Table S30a
estout B1 C1 D1 E1, style(tex) replace varwidth(25) collabels("") cells(b(star fmt(%9.2f)) se(par fmt(%9.2f))  p(par fmt(%9.8f))) stats(N, fmt(%9.2f %9.0f) labels("N")) starlevels(* 0.05 ** 0.01 *** 0.005) 

reg zika_neuro correction##ib4.response_quartile_exp if mosquito==0 & behavior==0, robust
est store B1

reg zika_bite correction##ib4.response_quartile_exp if mosquito==0 & behavior==0, robust
est store C1

reg zika_sex correction##ib4.response_quartile_exp if mosquito==0 & behavior==0, robust
est store D1

reg zika_handshake_belief correction##ib4.response_quartile_exp if mosquito==0 & behavior==0, robust
est store E1

*Table S31a
estout B1 C1 D1 E1, style(tex) replace varwidth(25) collabels("") cells(b(star fmt(%9.2f)) se(par fmt(%9.2f))  p(par fmt(%9.8f))) stats(N, fmt(%9.2f %9.0f) labels("N")) starlevels(* 0.05 ** 0.01 *** 0.005) 



/* STUDY 2 */

clear
set more off

use "zika-study2-yf.dta", clear

/* branch by disease */

count
destring random, replace
gen zika=(random==1)
gen yellow_fever=(random==2)
drop if random==.
count

/*descriptives*/

*gender
bysort zika: tab Q5 

*race
bysort zika: tab Q7

*age
bysort zika: tab Q4

*education
bysort zika: tab Q6

gen hs_degree=.
replace hs_degree=(Q6>18 & Q6!=20) if Q6!=.
bysort zika: tab hs_degree

*income
rename Q8 income
replace income=. if income==-99
bysort zika: su income, detail

* 1 is very left-leaning, 10 is very right-leaning *
rename Q16_1 ideology
bysort zika: su ideology, detail
bysort zika: tab ideology

gen left=(ideology<6) if ideology!=.
bysort zika: tab left

gen right=(ideology>5) if ideology!=.
bysort zika: tab right

recode Q125 (9=1 "Not accurate") (10=2 "A little accurate") (11=3 "More or less accurate") (12=4 "Very accurate") (-99=.), gen(cabals)
/* Major events like wars/recession/election are controlled by small groups of
people who are working secretly against the rest of us / the population */
/*
3	not recorded			-99
439	not accurate			9
536	a little accurate		10
941	more or less accurate	11
539	very accurate			12
*/

bysort zika: su cabals

/* Recode variables */

/* ZIKA control article */
gen ZikaControl=.
replace ZikaControl=0 if FL_18_DO_Condition2_Zikamyths_an==1 
replace ZikaControl=1 if FL_18_DO_Condition1_Zikacontrol==1

tab ZikaControl

/* ZIKA myths and facts article */
gen ZikaCorrection=.
replace ZikaCorrection=0 if FL_18_DO_Condition1_Zikacontrol==1
replace ZikaCorrection=1 if FL_18_DO_Condition2_Zikamyths_an==1 

/* YF control article */
gen YFControl=.
replace YFControl=0 if FL_56_DO_Condition4_YellowFeverm==1
replace YFControl=1 if FL_56_DO_Condition3_YellowFeverc==1

/* YF myths and facts article */
gen YFCorrection=.
replace YFCorrection=0 if FL_56_DO_Condition3_YellowFeverc==1
replace YFCorrection=1 if FL_56_DO_Condition4_YellowFeverm==1

egen sum = rowtotal(ZikaControl ZikaCorrection YFControl YFCorrection)
drop if sum==0

duplicates report psid

/* drop duplicates of id, keeping first finished observation in each group */

* tag duplicates *
duplicates tag psid, gen(duplicate)
list psid StartDate if duplicate>0
codebook psid
tab duplicate
duplicates list psid

* keep only the first duplicate *
bysort psid: egen mintime = min(StartDate)
format mintime %tc
gen keeper=(abs((StartDate-mintime)/1000)<=90)
sort psid StartDate
list psid StartDate mintime keeper Q4 Q65 Q204 if duplicate>0
drop if keeper==0 & duplicate>0
duplicates tag psid,gen(dupey2)
list psid StartDate mintime keeper Q4 Q65 Q204 if dupey2>0

tab ZikaControl ZikaCorrection if zika==1, missing
tab YFControl YFCorrection if zika==0, missing

/* response times*/

gen response_time=.
replace response_time=Q35_Page_Submit if ZikaCorrection==0 & zika==1
replace response_time=Q41_Page_Submit if ZikaCorrection==1 & zika==1
xtile response_quartile_exp=response_time,nq(4)

egen prior_response_time=rowtotal(Q108_Page_Submit Q225_Page_Submit Q226_Page_Submit Q106_Page_Submit Q109_Page_Submit Q111_Page_Submit Q118_Page_Submit Q119_Page_Submit Q112_Page_Submit Q114_Page_Submit Q116_Page_Submit Q117_Page_Submit)
xtile response_quartile=prior_response_time,nq(4)

/*trust*/

rename Q136 confidence_court 
/* How much confidence do you have in the Supreme Federal Court? */
/*
19	not recorded	-99
470	No confidence	1
428					2
428					3
489					4
418					5
208					6
110	Full confidence	7
*/
rename Q137 confidence_congress
/* How much confidence do you have in the National Congress? */
/*
16		not recorded	-99
1221	No confidence	1
540						2
347						3
255						4
113						5
43						6
35		Full confidence	7
*/
rename Q138 confidence_press
/* How much confidence do you have in the press? */
/*
11		not recorded	-99
365		No confidence	1
434						2
524						3
578						4
404						5
156						6
98		Full confidence	7
*/
rename Q139 confidence_govt
/* How much confidence do you have in the government? */
/*
20		not recorded	-99
1199	No confidence	1
580						2
338						3
264						4
107						5
31						6
31		Full confidence	7
*/
rename Q140 confidence_medicine
/* How much confidence do you have in medicines [doctors? portuguese word is "remedios"] ? */
/*
17		not recorded	-99
44		No confidence	1
116						2
277						3
549						4
749						5
554						6
264		Full confidence	7
*/
rename Q141 confidence_scientists
/* How much confidence do you have in scientists? */
/*
19		not recorded	-99
40		No confidence	1
102						2
168						3
396						4
718						5
705						6
422		Full confidence	7
*/
rename Q142 confidence_banks
/* How much confidence do you have in banks and financial institutions? */
/*
20		not recorded	-99
309		No confidence	1
400						2
489						3
604						4
459						5
209						6
80		Full confidence	7
*/
rename Q143 confidence_govhealth
/* How much confidence do you have in the people who manage the Ministry of 
Health */
/*
12		not recorded	-99
602		No confidence	1
564						2
514						3
474						4
260						5
95						6
49		Full confidence	7
*/

foreach var of varlist confidence* {
replace `var'=. if `var'==-99
}


/* Zika knowledge/CT DVs */

*Q45 - casual contact
recode Q45 (1=4) (2=3) (3=2) (4=1) (-99=.), gen(zika_handshake)
tab zika_handshake

*Q46 - GMOs cause Zika
recode Q46 (10=4) (11=3) (12=2) (13=1) (-99=.), gen(zika_gmo)
tab zika_gmo

*Q43 - Neurological problems
recode Q43 (6=1) (10=2) (11=3) (12=4) (-99=.), gen(zika_neuro)
tab zika_neuro

*Q47 - sexual transmission
recode Q47 (10=1) (11=2) (12=3) (13=4) (-99=.), gen(zika_sex)
tab zika_sex

*Q44 - transmitted by mosquitoes
recode Q44 (10=1) (11=2) (12=3) (13=4) (-99=.), gen(zika_bite)
tab zika_bite

*Q 148_NEW - weak immune systems
recode Q148_NEW (10=4) (11=3) (12=2) (13=1) (-99=.), gen(zika_immune)
tab zika_immune

*Q219 - larvicides cause microcephaly
recode Q219 (10=4) (11=3) (12=2) (13=1) (-99=.), gen(micro_larvicide)
tab micro_larvicide

*Q150_NEW - microcephaly increased b/c of gov’t mandating vaccines
recode Q150_NEW (10=4) (11=3) (12=2) (13=1) (-99=.), gen(micro_vaccine)
tab micro_vaccine

*Q151_NEW - transmission during pregnancy
recode Q151_NEW (10=1) (11=2) (12=3) (13=4) (-99=.), gen(zika_pregnancy)
tab zika_pregnancy

*Q152 - mother/baby microcephaly
recode Q152 (10=1) (11=2) (12=3) (13=4) (-99=.), gen(mom_zika)
tab mom_zika

*Q170 - transmission by blood 
recode Q170 (10=1) (11=2) (12=3) (13=4) (-99=.), gen(zika_blood)
tab zika_blood

*Q171 - microcephaly causes paralysis
recode Q171 (10=4) (11=3) (12=2) (13=1) (-99=.), gen(zika_paralysis)
tab zika_paralysis

/* Zika policy preference DVs */

*Q133 - larvicides

recode Q133 (4=1) (5=2) (6=3) (7=4) (8=5) (9=6) (10=7) (11=8) (12=9) (13=10) (-99=.), gen(policy_larvicide)
tab policy_larvicide

*Q53 - DTPa vax

recode Q53 (4=1) (5=2) (6=3) (7=4) (8=5) (9=6) (10=7) (11=8) (12=9) (13=10) (-99=.), gen(policy_dtap)
tab policy_dtap

*Q134 - use of GM mosquitoes
recode Q134 (4=1) (5=2) (6=3) (7=4) (8=5) (9=6) (10=7) (11=8) (12=9) (13=10) (-99=.), gen(policy_gmo)
tab policy_gmo

*Q135 - search of homes
recode Q135 (4=1) (5=2) (6=3) (7=4) (8=5) (9=6) (10=7) (11=8) (12=9) (13=10) (-99=.), gen(policy_homesearch)
tab policy_homesearch

/* Zika misc. DVs */

*Q64 - risk to human health
recode Q64 (7=1) (8=2) (9=3) (10=4) (12=5) (13=6) (-99=.), gen(risk_human)
tab risk_human

*Q65 - risk to personal health
recode Q65 (7=1) (8=2) (9=3) (10=4) (11=5) (12=6) (-99=.), gen(risk_personal)
tab risk_personal

*Q66 - concern about Zika
recode Q66 (1=1) (2=2) (3=3) (4=4) (5=5) (-99=.), gen(concern)
tab concern

*Q128 - wear long sleeves/pants
recode Q128 (4=5) (5=4) (6=3) (7=2) (8=1) (-99=.), gen(protect_sleeves)
tab protect_sleeves

*Q129 - use repellent
recode Q129 (4=5) (5=4) (6=3) (7=2) (8=1) (-99=.), gen(protect_repel)
tab protect_repel

*Q130 - screens/close windows
recode Q130 (4=5) (5=4) (6=3) (7=2) (8=1) (-99=.), gen(protect_screen)
tab protect_screen

/* Epistemic efficacy DVs */
recode Q221_1 (1=4) (2=3) (3=2) (4=1) (-99=.), gen(epe1_zika)
tab epe1_zika
recode Q221_4 (1=4) (2=3) (3=2) (4=1) (-99=.), gen(epe2_zika)
tab epe2_zika
recode Q221_5 (1=4) (2=3) (3=2) (4=1) (-99=.), gen(epe3_zika)
tab epe3_zika

/* Infection */

*Q69 - contracted Zika
recode Q69 (1=1) (2=0) (-99=.), gen(has_zika)
tab has_zika

*Q70 - friends/family contracted Zika
recode Q70 (1=1) (2=0) (-99=.), gen(family_zika)
tab family_zika

*Q215 - vaccinated
recode Q215 (1=1) (2=0) (-99=.), gen(vaccinated)
tab vaccinated

*Q217 - likelihood of getting vax in future
recode Q217 (1=1) (2=2) (3=3) (4=4) (5=5) (-99=.), gen(likely_vax)
tab likely_vax

/* YF knowledge/CT DVs */

*Q173 - transmitted by mosquitoes
recode Q173 (4=4) (3=3) (2=2) (1=1) (-99=.), gen(YF_bite)
tab YF_bite

*Q174 - transmitted by AA mosquito (same as Zika)
recode Q174 (4=4) (3=3) (2=2) (1=1) (-99=.), gen(YF_same)
tab YF_same

*Q175 - YF symptoms
recode Q175 (1=1) (2=2) (3=3) (4=4) (-99=.), gen(YF_symptoms)
tab YF_symptoms

*Q176 - YF can be fatal
recode Q176 (4=4) (3=3) (2=2) (1=1) (-99=.), gen(YF_fatal)
tab YF_fatal

*Q177 - country-wide vax recommendation
recode Q177 (4=4) (3=3) (2=2) (1=1) (-99=.), gen(YF_recommend)
tab YF_recommend

*Q178 - YF now common in cities
recode Q178 (4=4) (3=3) (2=2) (1=1) (-99=.), gen(YF_cities)
tab YF_cities

*Q179 - no effective YF vax
recode Q179 (1=4) (2=3) (3=2) (4=1) (-99=.), gen(YF_no_vax)
tab YF_no_vax

*Q180 - YF vax damages immune system
recode Q180 (1=4) (2=3) (3=2) (4=1) (-99=.), gen(YF_damage)
tab YF_damage

*Q181 - vax no longer effective b/c of genetic mutations
recode Q181 (1=4) (2=3) (3=2) (4=1) (-99=.), gen(YF_vax_ineffective)
tab YF_vax_ineffective

*Q182 - vax a fraud by big pharma
recode Q182 (1=4) (2=3) (3=2) (4=1) (-99=.), gen(YF_fraud)
tab YF_fraud

*Q183 - vax has life-threatening side effects
recode Q183 (1=4) (2=3) (3=2) (4=1) (-99=.), gen(YF_vax_effects)
tab YF_vax_effects

*Q184 - Propolis repeals mosquitoes
recode Q184 (1=4) (2=3) (3=2) (4=1) (-99=.), gen(YF_propolis)
tab YF_propolis

/* YF policy preference DVs */

*Q187 - YF vaccine to attend public schools
recode Q187 (1=1) (2=2) (3=3) (4=4) (5=5) (6=6) (7=7) (8=8) (9=9) (10=10) (-99=.), gen(YF_schools)
tab YF_schools

*Q188 - fines for non-vaccinators
recode Q188 (4=1) (5=2) (6=3) (7=4) (8=5) (9=6) (10=7) (11=8) (12=9) (13=10) (-99=.), gen(YF_fines)
tab YF_fines

*Q189 - use of larvicides in water
recode Q189 (4=1) (5=2) (6=3) (7=4) (8=5) (9=6) (10=7) (11=8) (12=9) (13=10) (-99=.), gen(YF_policy_larvicides)
tab YF_policy_larvicides

*Q190 - use of GMO mosquitoes
recode Q190 (4=1) (5=2) (6=3) (7=4) (8=5) (9=6) (10=7) (11=8) (12=9) (13=10) (-99=.), gen(YF_policy_gmo)
tab YF_policy_gmo

*Q191 - search of homes
recode Q191 (4=1) (5=2) (6=3) (7=4) (8=5) (9=6) (10=7) (11=8) (12=9) (13=10) (-99=.), gen(YF_policy_homesearch)
tab YF_policy_homesearch

/* YF misc. DVs */

*Q193 - YF risk to human health
recode Q193 (7=1) (8=2) (9=3) (10=4) (12=5) (13=6) (-99=.), gen(YF_risk_human)
tab YF_risk_human

*Q197 - YF risk to personal health
recode Q197 (7=1) (8=2) (9=3) (10=4) (11=5) (12=6) (-99=.), gen(YF_risk_personal)
tab YF_risk_personal

*Q199 - concern about YF
recode Q199 (1=1) (2=2) (3=3) (4=4) (5=5) (-99=.), gen(YF_concern)
tab YF_concern

*Q202 - wear long sleeves/pants
recode Q202 (4=5) (5=4) (6=3) (7=2) (8=1) (-99=.), gen(YF_protect_sleeves)
tab YF_protect_sleeves

*Q203 - use repellent
recode Q203 (4=5) (5=4) (6=3) (7=2) (8=1) (-99=.), gen(YF_protect_repel)
tab YF_protect_repel

*Q204 - screens/close windows
recode Q204 (4=5) (5=4) (6=3) (7=2) (8=1) (-99=.), gen(YF_protect_screen)
tab YF_protect_screen

/* Epistemic efficacy DVs */ 
recode Q223_1 (1=4) (2=3) (3=2) (4=1) (-99=.), gen(YF_epe1)
tab YF_epe1
recode Q223_4 (1=4) (2=3) (3=2) (4=1) (-99=.), gen(YF_epe2)
tab YF_epe2
recode Q223_5 (1=4) (2=3) (3=2) (4=1) (-99=.), gen(YF_epe3)
tab YF_epe3

/* Infection */

*Q206 - contracted Zika
recode Q206 (1=1) (2=0) (-99=.), gen(has_YF)
tab has_YF

*Q208 - friends/family contracted Zika
recode Q208 (1=1) (2=0) (-99=.), gen(family_YF)
tab family_YF

*Q210 - vaccinated
recode Q210 (1=1) (4=0) (-99=.), gen(YF_vaccinated)
tab YF_vaccinated

*Q212 - likelihood of getting vax in future
recode Q212 (1=1) (4=2) (5=3) (6=4) (7=5) (-99=.), gen(YF_likely_vax)
tab YF_likely_vax

/** Creating scales needed in analysis **/

*Zika knowledge 
factor zika_bite zika_neuro zika_sex zika_handshake zika_immune zika_pregnancy zika_blood zika_paralysis mom_zika, pcf factors(1) 
alpha zika_bite zika_neuro zika_sex zika_handshake zika_immune zika_pregnancy zika_blood zika_paralysis mom_zika, item casewise gen(ZikaKnow) asis

*Zika CTs 
factor zika_gmo micro_larvicide micro_vaccine, pcf factors(1) 
alpha zika_gmo micro_larvicide micro_vaccine, item casewise gen(ZikaMisp) 

*Zika policy prefs:
factor policy_gmo policy_dtap policy_larvicide policy_homesearch, pcf factors(1) 
alpha policy_gmo policy_dtap policy_larvicide policy_homesearch, item casewise gen(ZikaOpinions)

*Zika behavioral intentions:
factor protect_sleeves protect_repel protect_screen, pcf factors(1) 
alpha protect_sleeves protect_repel protect_screen, item casewise gen(ZikaBehavior) 

*YF knowledge:
factor YF_bite YF_same YF_symptoms YF_fatal YF_recommend YF_cities YF_no_vax YF_damage YF_fraud, pcf factors(1) 
alpha YF_bite YF_same YF_symptoms YF_fatal YF_recommend YF_cities YF_no_vax YF_damage YF_fraud, item casewise gen(YFKnow) 

*YF CTs:
factor YF_vax_ineffective YF_vax_effects YF_propolis, pcf factors(1) 
alpha YF_vax_ineffective YF_vax_effects YF_propolis, item casewise gen(YFMisp)

factor YF_schools YF_fines YF_policy_larvicides YF_policy_gmo YF_policy_homesearch, pcf factors(1) 
alpha YF_schools YF_fines YF_policy_larvicides YF_policy_gmo YF_policy_homesearch, item casewise gen(YFOpinions)

*YF behavioral intentions:
factor YF_protect_sleeves YF_protect_repel YF_protect_screen, pcf factors(1)
alpha YF_protect_sleeves YF_protect_repel YF_protect_screen, item casewise gen(YFBehavior)

*EPE:

forval i=1/3 {
gen epe`i'=.
replace epe`i'=epe`i'_zika if epe`i'_zika!=.
replace epe`i'=YF_epe`i' if YF_epe`i'!=.
}

su epe1-epe3

factor epe1 epe2 epe3, pcf factors(1)
alpha epe1 epe2 epe3, item casewise gen(EPE) 

/* Alternate coding of belief outcomes (higher=more belief) */

/* Zika */

*Q45 - casual contact
recode Q45 (1=1) (2=2) (3=3) (4=4) (-99=.), gen(zika_handshake_belief)
tab zika_handshake_belief

*Q46 - GMOs cause Zika
recode Q46 (10=1) (11=2) (12=3) (13=4) (-99=.), gen(zika_gmo_belief)
tab zika_gmo_belief

*Q43 - Neurological problems
recode Q43 (6=1) (10=2) (11=3) (12=4) (-99=.), gen(zika_neuro_belief)
tab zika_neuro_belief

*Q47 - sexual transmission
recode Q47 (10=1) (11=2) (12=3) (13=4) (-99=.), gen(zika_sex_belief)
tab zika_sex_belief

*Q44 - transmitted by mosquitoes
recode Q44 (10=1) (11=2) (12=3) (13=4) (-99=.), gen(zika_bite_belief)
tab zika_bite_belief

*Q148_NEW - weak immune systems
recode Q148_NEW (10=1) (11=2) (12=3) (13=4) (-99=.), gen(zika_immune_belief)
tab zika_immune_belief

*Q219 - larvicides cause microcephaly
recode Q219 (10=1) (11=2) (12=3) (13=4) (-99=.), gen(micro_larvicide_belief)
tab micro_larvicide_belief

*Q150_NEW - microcephaly increased b/c of gov’t mandating vaccines
recode Q150_NEW (10=1) (11=2) (12=3) (13=4) (-99=.), gen(micro_vaccine_belief)
tab micro_vaccine_belief

*Q151_NEW - transmission during pregnancy
recode Q151_NEW (10=1) (11=2) (12=3) (13=4) (-99=.), gen(zika_pregnancy_belief)
tab zika_pregnancy_belief

*Q152 - mother/baby microcephaly
recode Q152 (10=1) (11=2) (12=3) (13=4) (-99=.), gen(mom_zika_belief)
tab mom_zika_belief

*Q170 - transmission by blood 
recode Q170 (10=1) (11=2) (12=3) (13=4) (-99=.), gen(zika_blood_belief)
tab zika_blood_belief

*Q171 - microcephaly causes paralysis
recode Q171 (10=1) (11=2) (12=3) (13=4) (-99=.), gen(zika_paralysis_belief)
tab zika_paralysis_belief

/* YF */

*Q173 - transmitted by mosquitoes
recode Q173 (4=4) (3=3) (2=2) (1=1) (-99=.), gen(YF_bite_belief)
tab YF_bite_belief

*Q174 - transmitted by AA mosquito (same as Zika)
recode Q174 (4=4) (3=3) (2=2) (1=1) (-99=.), gen(YF_same_belief)
tab YF_same_belief

*Q175 - YF symptoms
recode Q175 (1=1) (2=2) (3=3) (4=4) (-99=.), gen(YF_symptoms_belief)
tab YF_symptoms_belief

*Q176 - YF can be fatal
recode Q176 (4=4) (3=3) (2=2) (1=1) (-99=.), gen(YF_fatal_belief)
tab YF_fatal_belief

*Q177 - country-wide vax recommendation
recode Q177 (4=4) (3=3) (2=2) (1=1) (-99=.), gen(YF_recommend_belief)
tab YF_recommend_belief

*Q178 - YF now common in cities
recode Q178 (4=4) (3=3) (2=2) (1=1) (-99=.), gen(YF_cities_belief)
tab YF_cities_belief

*Q179 - no effective YF vax
recode Q179 (1=1) (2=2) (3=3) (4=4) (-99=.), gen(YF_no_vax_belief)
tab YF_no_vax_belief

*Q180 - YF vax damages immune system
recode Q180 (1=1) (2=2) (3=3) (4=4) (-99=.), gen(YF_damage_belief)
tab YF_damage_belief

*Q181 - vax no longer effective b/c of genetic mutations
recode Q181 (1=1) (2=2) (3=3) (4=4) (-99=.), gen(YF_vax_ineffective_belief)
tab YF_vax_ineffective_belief

*Q182 - vax a fraud by big pharma
recode Q182 (1=1) (2=2) (3=3) (4=4) (-99=.), gen(YF_fraud_belief)
tab YF_fraud_belief

*Q183 - vax has life-threatening side effects
recode Q183 (1=1) (2=2) (3=3) (4=4) (-99=.), gen(YF_vax_effects_belief)
tab YF_vax_effects_belief

*Q184 - Propolis repeals mosquitoes
recode Q184 (1=1) (2=2) (3=3) (4=4) (-99=.), gen(YF_propolis_belief)
tab YF_propolis_belief

/** Creating scales of BELIEF CODED OUTCOMES needed in analysis of H5 **/

*Zika knowledge belief
factor zika_bite_belief zika_neuro_belief zika_sex_belief zika_handshake_belief zika_immune_belief zika_pregnancy_belief zika_blood_belief zika_paralysis_belief mom_zika_belief, pcf factors(1) 
alpha zika_bite_belief zika_neuro_belief zika_sex_belief zika_handshake_belief zika_immune_belief zika_pregnancy_belief zika_blood_belief zika_paralysis_belief mom_zika_belief, item casewise gen(ZikaKnow_belief) asis  

alpha zika_bite_belief zika_neuro_belief zika_sex_belief zika_pregnancy_belief zika_blood_belief mom_zika_belief, item casewise gen(ZikaKnow_belief_true) asis  

*YF knowledge belief

factor YF_bite_belief YF_same_belief YF_symptoms_belief YF_fatal_belief YF_recommend_belief YF_cities_belief YF_no_vax_belief YF_damage_belief YF_fraud_belief, pcf factors(1) 

alpha YF_bite_belief YF_same_belief YF_symptoms_belief YF_fatal_belief YF_recommend_belief YF_cities_belief YF_no_vax_belief YF_damage_belief YF_fraud_belief, item casewise gen(YFKnow_belief) asis 

alpha YF_bite_belief YF_same_belief YF_symptoms_belief YF_fatal_belief YF_recommend_belief YF_cities_belief, item casewise gen(YFKnow_belief_true) asis 

*Zika CT belief

factor zika_gmo_belief micro_larvicide_belief micro_vaccine_belief, pcf factors(1) 
alpha zika_gmo_belief micro_larvicide_belief micro_vaccine_belief, item casewise gen(ZikaMisp_belief) 

*YF CT belief

factor YF_vax_ineffective_belief YF_vax_effects_belief YF_propolis_belief, pcf factors(1) 
alpha YF_vax_ineffective_belief YF_vax_effects_belief YF_propolis_belief, item casewise gen(YFMisp_belief) 

save "ZikaYF_study 2-cleaned.dta", replace


/* Analysis */

reg ZikaMisp_belief ZikaCorrection, robust
est store A

reg zika_gmo_belief ZikaCorrection, robust
est store B

reg micro_larvicide_belief ZikaCorrection, robust
est store C

reg micro_vaccine_belief ZikaCorrection, robust
est store D

*Table 2b
estout A B C D, style(tex) replace varwidth(25) collabels("") cells(b(star fmt(%9.2f)) se(par fmt(%9.2f))) stats(N, fmt(%9.2f %9.0f) labels("N")) starlevels(* 0.05 ** 0.01 *** 0.005) 

oprobit zika_gmo_belief ZikaCorrection, robust
est store B

oprobit micro_larvicide_belief ZikaCorrection, robust
est store C

oprobit micro_vaccine_belief ZikaCorrection, robust
est store D

*Table S14a
estout B C D, style(tex) replace varwidth(25) collabels("") cells(b(star fmt(%9.2f)) se(par fmt(%9.2f))) stats(N, fmt(%9.2f %9.0f) labels("N")) starlevels(* 0.05 ** 0.01 *** 0.005) 

/*related*/
reg zika_bite_belief ZikaCorrection, robust
est store B

/*unrelated (h5a)*/
reg zika_neuro_belief ZikaCorrection, robust
est store A

reg zika_sex_belief ZikaCorrection, robust
est store C

reg zika_handshake_belief ZikaCorrection, robust
est store D

reg zika_immune_belief ZikaCorrection, robust
est store E

reg zika_pregnancy_belief ZikaCorrection, robust
est store F

reg mom_zika_belief ZikaCorrection, robust
est store G

reg zika_blood_belief ZikaCorrection, robust
est store H

reg zika_paralysis_belief ZikaCorrection, robust
est store I

*Table 3b
estout A B C D E F G H I, style(tex) replace varwidth(25) collabels("") cells(b(star fmt(%9.2f)) se(par fmt(%9.2f))) stats(N, fmt(%9.2f %9.0f) labels("N")) starlevels(* 0.05 ** 0.01 *** 0.005) 

*joint test
mvreg zika_neuro_belief zika_bite_belief zika_sex_belief zika_handshake_belief zika_immune_belief zika_pregnancy_belief mom_zika_belief zika_blood_belief zika_paralysis_belief = ZikaCorrection
test ZikaCorrection

*effect sizes

foreach var of varlist zika_bite_belief zika_neuro_belief zika_sex_belief zika_handshake_belief zika_immune_belief zika_pregnancy_belief zika_blood_belief zika_paralysis_belief mom_zika_belief {
esize twosample `var' if ZikaControl==1 | ZikaCorrection==1, by(ZikaCorrection) cohensd 
reg `var' ZikaCorrection, robust
}

oprobit zika_bite_belief ZikaCorrection, robust
est store B

oprobit zika_neuro_belief ZikaCorrection, robust
est store A

oprobit zika_sex_belief ZikaCorrection, robust
est store C

oprobit zika_handshake_belief ZikaCorrection, robust
est store D

oprobit zika_immune_belief ZikaCorrection, robust
est store E

oprobit zika_pregnancy_belief ZikaCorrection, robust
est store F

oprobit mom_zika_belief ZikaCorrection, robust
est store G

oprobit zika_blood_belief ZikaCorrection, robust
est store H

oprobit zika_paralysis_belief ZikaCorrection, robust
est store I

*Table S14b
estout A B C D E F G H I, style(tex) replace varwidth(25) collabels("") cells(b(star fmt(%9.2f)) se(par fmt(%9.2f))) stats(N, fmt(%9.2f %9.0f) labels("N")) starlevels(* 0.05 ** 0.01 *** 0.005) 

reg YFMisp_belief YFCorrection, robust
est store A 

reg YF_vax_ineffective_belief YFCorrection, robust
est store B

reg YF_vax_effects_belief YFCorrection, robust
est store C

reg YF_propolis_belief YFCorrection, robust
est store D

*Table 4a
estout A B C D, style(tex) replace varwidth(25) collabels("") cells(b(star fmt(%9.2f)) se(par fmt(%9.2f))) stats(N, fmt(%9.2f %9.0f) labels("N")) starlevels(* 0.05 ** 0.01 *** 0.005) 

oprobit YF_vax_ineffective_belief YFCorrection, robust
est store B

oprobit YF_vax_effects_belief YFCorrection, robust
est store C

oprobit YF_propolis_belief YFCorrection, robust
est store D

*Table S23a
estout B C D, style(tex) replace varwidth(25) collabels("") cells(b(star fmt(%9.2f)) se(par fmt(%9.2f))) stats(N, fmt(%9.2f %9.0f) labels("N")) starlevels(* 0.05 ** 0.01 *** 0.005) 

reg YF_bite_belief YFCorrection, robust
est store A

reg YF_no_vax_belief YFCorrection, robust
est store B

reg YF_same_belief YFCorrection, robust
est store C

reg YF_symptoms_belief YFCorrection, robust
est store D

reg YF_fatal_belief YFCorrection, robust
est store E

reg YF_recommend_belief YFCorrection, robust
est store F

reg YF_cities_belief YFCorrection, robust
est store G

reg YF_damage_belief YFCorrection, robust
est store H

reg YF_fraud_belief YFCorrection, robust
est store I

*Table 4b
estout A B C D E F G H I, style(tex) replace varwidth(25) collabels("") cells(b(star fmt(%9.2f)) se(par fmt(%9.2f))) stats(N, fmt(%9.2f %9.0f) labels("N")) starlevels(* 0.05 ** 0.01 *** 0.005) 

oprobit YF_bite_belief YFCorrection, robust
est store A

oprobit YF_no_vax_belief YFCorrection, robust
est store B

oprobit YF_same_belief YFCorrection, robust 
est store C

oprobit YF_symptoms_belief YFCorrection, robust
est store D

oprobit YF_fatal_belief YFCorrection, robust
est store E

oprobit YF_recommend_belief YFCorrection, robust
est store F

oprobit YF_cities_belief YFCorrection, robust
est store G

oprobit YF_damage_belief YFCorrection, robust
est store H

oprobit YF_fraud_belief YFCorrection, robust
est store I

*Table S23b
estout A B C D E F G H I, style(tex) replace varwidth(25) collabels("") cells(b(star fmt(%9.2f)) se(par fmt(%9.2f))) stats(N, fmt(%9.2f %9.0f) labels("N")) starlevels(* 0.05 ** 0.01 *** 0.005) 

reg ZikaOpinions ZikaCorrection, robust
est store A

reg ZikaBehavior ZikaCorrection, robust
est store B

*Table S13
estout A B, style(tex) replace varwidth(25) collabels("") cells(b(star fmt(%9.2f)) se(par fmt(%9.2f))) stats(N, fmt(%9.2f %9.0f) labels("N")) starlevels(* 0.05 ** 0.01 *** 0.005) 

reg YFOpinions YFCorrection, robust 
est store A1

reg YFBehavior YFCorrection, robust
est store B1

*Table S22
estout A1 B1, style(tex) replace varwidth(25) collabels("") cells(b(star fmt(%9.2f)) se(par fmt(%9.2f))) stats(N, fmt(%9.2f %9.0f) labels("N")) starlevels(* 0.05 ** 0.01 *** 0.005) 

gen college=0
replace college=1 if Q6==23

reg zika_bite_belief c.college##ZikaCorrection if ZikaControl==1 | ZikaCorrection==1, robust
est store B

reg zika_neuro_belief c.college##ZikaCorrection if ZikaControl==1 | ZikaCorrection==1, robust
est store A

reg zika_sex_belief c.college##ZikaCorrection if ZikaControl==1 | ZikaCorrection==1, robust
est store C

reg zika_handshake_belief c.college##ZikaCorrection if ZikaControl==1 | ZikaCorrection==1, robust
est store D

reg zika_immune_belief c.college##ZikaCorrection if ZikaControl==1 | ZikaCorrection==1, robust
est store E

reg zika_pregnancy_belief c.college##ZikaCorrection if ZikaControl==1 | ZikaCorrection==1, robust
est store F

reg mom_zika_belief c.college##ZikaCorrection if ZikaControl==1 | ZikaCorrection==1, robust
est store G

reg zika_blood_belief c.college##ZikaCorrection if ZikaControl==1 | ZikaCorrection==1, robust
est store H

reg zika_paralysis_belief c.college##ZikaCorrection if ZikaControl==1 | ZikaCorrection==1, robust
est store I

*Table S20a
estout A B C D E F G H I, style(tex) replace varwidth(25) collabels("") cells(b(star fmt(%9.2f)) se(par fmt(%9.2f))) stats(N, fmt(%9.2f %9.0f) labels("N")) starlevels(* 0.05 ** 0.01 *** 0.005) 

gen sciknow1=(Q23_1==1)
gen sciknow2=(Q23_2==2)
gen sciknow3=(Q23_3==2)
gen sciknow4=(Q23_4==1)
gen sciknow5=(Q23_5==1)
gen sciknow6=(Q23_6==2)
gen sciknow7=(Q25==1)
gen sciknow8=(Q24==3)

factor sciknow*, pcf
rotate, varimax

alpha sciknow*, item 
egen sciknow=rowtotal(sciknow*)
replace sciknow=sciknow/8

reg zika_bite_belief c.sciknow##ZikaCorrection if ZikaControl==1 | ZikaCorrection==1, robust
est store B

reg zika_neuro_belief c.sciknow##ZikaCorrection if ZikaControl==1 | ZikaCorrection==1, robust
est store A

reg zika_sex_belief c.sciknow##ZikaCorrection if ZikaControl==1 | ZikaCorrection==1, robust
est store C

reg zika_handshake_belief c.sciknow##ZikaCorrection if ZikaControl==1 | ZikaCorrection==1, robust
est store D

reg zika_immune_belief c.sciknow##ZikaCorrection if ZikaControl==1 | ZikaCorrection==1, robust
est store E

reg zika_pregnancy_belief c.sciknow##ZikaCorrection if ZikaControl==1 | ZikaCorrection==1, robust
est store F

reg mom_zika_belief c.sciknow##ZikaCorrection if ZikaControl==1 | ZikaCorrection==1, robust
est store G

reg zika_blood_belief c.sciknow##ZikaCorrection if ZikaControl==1 | ZikaCorrection==1, robust
est store H

reg zika_paralysis_belief c.sciknow##ZikaCorrection if ZikaControl==1 | ZikaCorrection==1, robust
est store I

*Table S20b
estout A B C D E F G H I, style(tex) replace varwidth(25) collabels("") cells(b(star fmt(%9.2f)) se(par fmt(%9.2f))) stats(N, fmt(%9.2f %9.0f) labels("N")) starlevels(* 0.05 ** 0.01 *** 0.005) 

reg zika_neuro_belief ZikaCorrection##c.cabal, robust
lincom 1.ZikaCorrection+2*1.ZikaCorrection#cabal
lincom 1.ZikaCorrection+3*1.ZikaCorrection#cabal
lincom 1.ZikaCorrection+4*1.ZikaCorrection#cabal
est store A1

reg zika_neuro_belief confidence_govt ZikaCorrection##c.confidence_govt, robust
est store A2

reg zika_neuro_belief ZikaCorrection##c.confidence_govhealth, robust
est store A3

reg zika_neuro_belief ZikaCorrection##c.confidence_medicine, robust
est store A4

reg zika_neuro_belief ZikaCorrection##c.confidence_scientists, robust
est store A5

reg zika_bite_belief ZikaCorrection##c.cabal, robust
lincom 1.ZikaCorrection+2*1.ZikaCorrection#cabal
lincom 1.ZikaCorrection+3*1.ZikaCorrection#cabal
lincom 1.ZikaCorrection+4*1.ZikaCorrection#cabal
est store B1

reg zika_bite_belief confidence_govt ZikaCorrection##c.confidence_govt, robust
est store B2

reg zika_bite_belief ZikaCorrection##c.confidence_govhealth, robust
est store B3

reg zika_bite_belief ZikaCorrection##c.confidence_medicine, robust
est store B4

reg zika_bite_belief ZikaCorrection##c.confidence_scientists, robust
est store B5

reg zika_handshake_belief ZikaCorrection##c.cabal, robust
lincom 1.ZikaCorrection+2*1.ZikaCorrection#cabal
lincom 1.ZikaCorrection+3*1.ZikaCorrection#cabal
lincom 1.ZikaCorrection+4*1.ZikaCorrection#cabal
est store D1

reg zika_handshake_belief confidence_govt ZikaCorrection##c.confidence_govt, robust
est store D2

reg zika_handshake_belief ZikaCorrection##c.confidence_govhealth, robust
est store D3

reg zika_handshake_belief ZikaCorrection##c.confidence_medicine, robust
est store D4

reg zika_handshake_belief ZikaCorrection##c.confidence_scientists, robust
est store D5

reg zika_sex_belief ZikaCorrection##c.cabal, robust
lincom 1.ZikaCorrection+2*1.ZikaCorrection#cabal
lincom 1.ZikaCorrection+3*1.ZikaCorrection#cabal
lincom 1.ZikaCorrection+4*1.ZikaCorrection#cabal
est store C1

reg zika_sex_belief confidence_govt ZikaCorrection##c.confidence_govt, robust
est store C2

reg zika_sex_belief ZikaCorrection##c.confidence_govhealth, robust
est store C3

reg zika_sex_belief ZikaCorrection##c.confidence_medicine, robust
est store C4

reg zika_sex_belief ZikaCorrection##c.confidence_scientists, robust
est store C5

reg zika_immune_belief ZikaCorrection##c.cabal, robust
lincom 1.ZikaCorrection+2*1.ZikaCorrection#cabal
lincom 1.ZikaCorrection+3*1.ZikaCorrection#cabal
lincom 1.ZikaCorrection+4*1.ZikaCorrection#cabal
est store A1b

reg zika_immune_belief confidence_govt ZikaCorrection##c.confidence_govt, robust
est store A1b2

reg zika_immune_belief ZikaCorrection##c.confidence_govhealth, robust
est store A1b3

reg zika_immune_belief ZikaCorrection##c.confidence_medicine, robust
est store A1b4

reg zika_immune_belief ZikaCorrection##c.confidence_scientists, robust
est store A1b5

reg zika_pregnancy_belief ZikaCorrection##c.cabal, robust
lincom 1.ZikaCorrection+2*1.ZikaCorrection#cabal
lincom 1.ZikaCorrection+3*1.ZikaCorrection#cabal
lincom 1.ZikaCorrection+4*1.ZikaCorrection#cabal
est store A1c

reg zika_pregnancy_belief confidence_govt ZikaCorrection##c.confidence_govt, robust
est store A1c2

reg zika_pregnancy_belief ZikaCorrection##c.confidence_govhealth, robust
est store A1c3

reg zika_pregnancy_belief ZikaCorrection##c.confidence_medicine, robust
est store A1c4

reg zika_pregnancy_belief ZikaCorrection##c.confidence_scientists, robust
est store A1c5

reg zika_blood_belief ZikaCorrection##c.cabal, robust
lincom 1.ZikaCorrection+2*1.ZikaCorrection#cabal
lincom 1.ZikaCorrection+3*1.ZikaCorrection#cabal
lincom 1.ZikaCorrection+4*1.ZikaCorrection#cabal
est store A1e 

reg zika_blood_belief confidence_govt ZikaCorrection##c.confidence_govt, robust
est store A1e2

reg zika_blood_belief ZikaCorrection##c.confidence_govhealth, robust
est store A1e3

reg zika_blood_belief ZikaCorrection##c.confidence_medicine, robust
est store A1e4

reg zika_blood_belief ZikaCorrection##c.confidence_scientists, robust
est store A1e5

reg zika_paralysis_belief ZikaCorrection##c.cabal, robust
lincom 1.ZikaCorrection+2*1.ZikaCorrection#cabal
lincom 1.ZikaCorrection+3*1.ZikaCorrection#cabal
lincom 1.ZikaCorrection+4*1.ZikaCorrection#cabal
est store A1f

reg zika_paralysis_belief confidence_govt ZikaCorrection##c.confidence_govt, robust
est store A1f2

reg zika_paralysis_belief ZikaCorrection##c.confidence_govhealth, robust
est store A1f3

reg zika_paralysis_belief ZikaCorrection##c.confidence_medicine, robust
est store A1f4

reg zika_paralysis_belief ZikaCorrection##c.confidence_scientists, robust
est store A1f5

reg mom_zika_belief ZikaCorrection##c.cabal, robust
lincom 1.ZikaCorrection+2*1.ZikaCorrection#cabal
lincom 1.ZikaCorrection+3*1.ZikaCorrection#cabal
lincom 1.ZikaCorrection+4*1.ZikaCorrection#cabal
est store A1d

reg mom_zika_belief confidence_govt ZikaCorrection##c.confidence_govt, robust
est store A1d2

reg mom_zika_belief ZikaCorrection##c.confidence_govhealth, robust
est store A1d3

reg mom_zika_belief ZikaCorrection##c.confidence_medicine, robust
est store A1d4

reg mom_zika_belief ZikaCorrection##c.confidence_scientists, robust
est store A1d5

reg ZikaMisp ZikaCorrection##c.cabal, robust
lincom 1.ZikaCorrection+2*1.ZikaCorrection#cabal
lincom 1.ZikaCorrection+3*1.ZikaCorrection#cabal
lincom 1.ZikaCorrection+4*1.ZikaCorrection#cabal
est store E1

reg ZikaMisp confidence_govt ZikaCorrection##c.confidence_govt, robust
est store E2

reg ZikaMisp ZikaCorrection##c.confidence_govhealth, robust
est store E3

reg ZikaMisp ZikaCorrection##c.confidence_medicine, robust
est store E4

reg ZikaMisp ZikaCorrection##c.confidence_scientists, robust
est store E5

reg ZikaOpinions ZikaCorrection##c.cabal, robust 
lincom 1.ZikaCorrection+2*1.ZikaCorrection#cabal
lincom 1.ZikaCorrection+3*1.ZikaCorrection#cabal
lincom 1.ZikaCorrection+4*1.ZikaCorrection#cabal
est store F1

reg ZikaOpinions ZikaCorrection##c.confidence_govt, robust
est store F2

reg ZikaOpinions ZikaCorrection##c.confidence_govhealth, robust
est store F3

reg ZikaOpinions ZikaCorrection##c.confidence_medicine, robust
est store F4

reg ZikaOpinions ZikaCorrection##c.confidence_scientists, robust
est store F5

reg ZikaBehavior ZikaCorrection##c.cabal, robust 
lincom 1.ZikaCorrection+2*1.ZikaCorrection#cabal
lincom 1.ZikaCorrection+3*1.ZikaCorrection#cabal
lincom 1.ZikaCorrection+4*1.ZikaCorrection#cabal
est store G1

reg ZikaBehavior ZikaCorrection##c.confidence_govt, robust
est store G2

reg ZikaBehavior ZikaCorrection##c.confidence_govhealth, robust
est store G3

reg ZikaBehavior ZikaCorrection##c.confidence_medicine, robust
est store G4

reg ZikaBehavior ZikaCorrection##c.confidence_scientists, robust
est store G5

*Table S15
estout A1 B1 C1 D1 A1b A1c A1d A1e A1f E1 F1 G1, style(tex) replace varwidth(25) collabels("") cells(b(star fmt(%9.2f)) se(par fmt(%9.2f))) stats(N, fmt(%9.2f %9.0f) labels("N")) starlevels(* 0.05 ** 0.01 *** 0.005) 

*Table S16
estout A2 B2 C2 D2 A1b2 A1c2 A1d2 A1e2 A1f2 E2 F2 G2, style(tex) replace varwidth(25) collabels("") cells(b(star fmt(%9.2f)) se(par fmt(%9.2f))) stats(N, fmt(%9.2f %9.0f) labels("N")) starlevels(* 0.05 ** 0.01 *** 0.005) 

*Table S17
estout A3 B3 C3 D3 A1b3 A1c3 A1d3 A1e3 A1f3 E3 F3 G3, style(tex) replace varwidth(25) collabels("") cells(b(star fmt(%9.2f)) se(par fmt(%9.2f))) stats(N, fmt(%9.2f %9.0f) labels("N")) starlevels(* 0.05 ** 0.01 *** 0.005) 

*Table S18
estout A4 B4 C4 D4 A1b4 A1c4 A1d4 A1e4 A1f4 E4 F4 G4, style(tex) replace varwidth(25) collabels("") cells(b(star fmt(%9.2f)) se(par fmt(%9.2f))) stats(N, fmt(%9.2f %9.0f) labels("N")) starlevels(* 0.05 ** 0.01 *** 0.005) 

*Table S19
estout A5 B5 C5 D5 A1b5 A1c5 A1d5 A1e5 A1f5 E5 F5 G5, style(tex) replace varwidth(25) collabels("") cells(b(star fmt(%9.2f)) se(par fmt(%9.2f))) stats(N, fmt(%9.2f %9.0f) labels("N")) starlevels(* 0.05 ** 0.01 *** 0.005) 

reg YF_bite_belief YFCorrection##c.cabal, robust
lincom 1.YFCorrection+2*1.YFCorrection#cabal
lincom 1.YFCorrection+3*1.YFCorrection#cabal
lincom 1.YFCorrection+4*1.YFCorrection#cabal
est store A1

reg YF_bite_belief confidence_govt YFCorrection##c.confidence_govt, robust
est store A2

reg YF_bite_belief YFCorrection##c.confidence_govhealth, robust
est store A3

reg YF_bite_belief YFCorrection##c.confidence_medicine, robust
est store A4

reg YF_bite_belief YFCorrection##c.confidence_scientists, robust
est store A5

reg YF_no_vax_belief YFCorrection##c.cabal, robust
lincom 1.YFCorrection+2*1.YFCorrection#cabal
lincom 1.YFCorrection+3*1.YFCorrection#cabal
lincom 1.YFCorrection+4*1.YFCorrection#cabal
est store B1

reg YF_no_vax_belief confidence_govt YFCorrection##c.confidence_govt, robust
est store B2

reg YF_no_vax_belief YFCorrection##c.confidence_govhealth, robust
est store B3

reg YF_no_vax_belief YFCorrection##c.confidence_medicine, robust
est store B4

reg YF_no_vax_belief YFCorrection##c.confidence_scientists, robust
est store B5

reg YF_symptoms_belief YFCorrection##c.cabal, robust
lincom 1.YFCorrection+2*1.YFCorrection#cabal
lincom 1.YFCorrection+3*1.YFCorrection#cabal
lincom 1.YFCorrection+4*1.YFCorrection#cabal
est store D1

reg YF_symptoms_belief confidence_govt YFCorrection##c.confidence_govt, robust
est store D2

reg YF_symptoms_belief YFCorrection##c.confidence_govhealth, robust
est store D3

reg YF_symptoms_belief YFCorrection##c.confidence_medicine, robust
est store D4

reg YF_symptoms_belief YFCorrection##c.confidence_scientists, robust
est store D5

reg YF_same_belief YFCorrection##c.cabal, robust
lincom 1.YFCorrection+2*1.YFCorrection#cabal
lincom 1.YFCorrection+3*1.YFCorrection#cabal
lincom 1.YFCorrection+4*1.YFCorrection#cabal
est store C1

reg YF_same_belief confidence_govt YFCorrection##c.confidence_govt, robust
est store C2

reg YF_same_belief YFCorrection##c.confidence_govhealth, robust
est store C3

reg YF_same_belief YFCorrection##c.confidence_medicine, robust
est store C4

reg YF_same_belief YFCorrection##c.confidence_scientists, robust
est store C5

reg YF_fatal_belief YFCorrection##c.cabal, robust
lincom 1.YFCorrection+2*1.YFCorrection#cabal
lincom 1.YFCorrection+3*1.YFCorrection#cabal
lincom 1.YFCorrection+4*1.YFCorrection#cabal
est store A1b

reg YF_fatal_belief confidence_govt YFCorrection##c.confidence_govt, robust
est store A1b2

reg YF_fatal_belief YFCorrection##c.confidence_govhealth, robust
est store A1b3

reg YF_fatal_belief YFCorrection##c.confidence_medicine, robust
est store A1b4

reg YF_fatal_belief YFCorrection##c.confidence_scientists, robust
est store A1b5

reg YF_recommend_belief YFCorrection##c.cabal, robust
lincom 1.YFCorrection+2*1.YFCorrection#cabal
lincom 1.YFCorrection+3*1.YFCorrection#cabal
lincom 1.YFCorrection+4*1.YFCorrection#cabal
est store A1c

reg YF_recommend_belief confidence_govt YFCorrection##c.confidence_govt, robust
est store A1c2

reg YF_recommend_belief YFCorrection##c.confidence_govhealth, robust
est store A1c3

reg YF_recommend_belief YFCorrection##c.confidence_medicine, robust
est store A1c4

reg YF_recommend_belief YFCorrection##c.confidence_scientists, robust
est store A1c5

reg YF_cities_belief YFCorrection##c.cabal, robust
lincom 1.YFCorrection+2*1.YFCorrection#cabal
lincom 1.YFCorrection+3*1.YFCorrection#cabal
lincom 1.YFCorrection+4*1.YFCorrection#cabal
est store A1d 

reg YF_cities_belief confidence_govt YFCorrection##c.confidence_govt, robust
est store A1d2

reg YF_cities_belief YFCorrection##c.confidence_govhealth, robust
est store A1d3

reg YF_cities_belief YFCorrection##c.confidence_medicine, robust
est store A1d4

reg YF_cities_belief YFCorrection##c.confidence_scientists, robust
est store A1d5

reg YF_damage_belief YFCorrection##c.cabal, robust
lincom 1.YFCorrection+2*1.YFCorrection#cabal
lincom 1.YFCorrection+3*1.YFCorrection#cabal
lincom 1.YFCorrection+4*1.YFCorrection#cabal
est store A1e

reg YF_damage_belief confidence_govt YFCorrection##c.confidence_govt, robust
est store A1e2

reg YF_damage_belief YFCorrection##c.confidence_govhealth, robust
est store A1e3

reg YF_damage_belief YFCorrection##c.confidence_medicine, robust
est store A1e4

reg YF_damage_belief YFCorrection##c.confidence_scientists, robust
est store A1e5

reg YF_fraud_belief YFCorrection##c.cabal, robust
lincom 1.YFCorrection+2*1.YFCorrection#cabal
lincom 1.YFCorrection+3*1.YFCorrection#cabal
lincom 1.YFCorrection+4*1.YFCorrection#cabal
est store A1f

reg YF_fraud_belief confidence_govt YFCorrection##c.confidence_govt, robust
est store A1f2

reg YF_fraud_belief YFCorrection##c.confidence_govhealth, robust
est store A1f3

reg YF_fraud_belief YFCorrection##c.confidence_medicine, robust
est store A1f4

reg YF_fraud_belief YFCorrection##c.confidence_scientists, robust
est store A1f5

reg YFMisp YFCorrection##c.cabal, robust
lincom 1.YFCorrection+2*1.YFCorrection#cabal
lincom 1.YFCorrection+3*1.YFCorrection#cabal
lincom 1.YFCorrection+4*1.YFCorrection#cabal
est store E1

reg YFMisp confidence_govt YFCorrection##c.confidence_govt, robust
est store E2

reg YFMisp YFCorrection##c.confidence_govhealth, robust
est store E3

reg YFMisp YFCorrection##c.confidence_medicine, robust
est store E4

reg YFMisp YFCorrection##c.confidence_scientists, robust
est store E5

reg YFOpinions cabal YFCorrection##c.cabal, robust 
lincom 1.YFCorrection+2*1.YFCorrection#cabal
lincom 1.YFCorrection+3*1.YFCorrection#cabal
lincom 1.YFCorrection+4*1.YFCorrection#cabal
est store F1

reg YFOpinions YFCorrection##c.confidence_govt, robust
est store F2

reg YFOpinions YFCorrection##c.confidence_govhealth, robust
est store F3

reg YFOpinions YFCorrection##c.confidence_medicine, robust
est store F4

reg YFOpinions YFCorrection##c.confidence_scientists, robust
est store F5

reg YFBehavior cabal YFCorrection##c.cabal, robust 
lincom 1.YFCorrection+2*1.YFCorrection#cabal
lincom 1.YFCorrection+3*1.YFCorrection#cabal
lincom 1.YFCorrection+4*1.YFCorrection#cabal
est store G1

reg YFBehavior YFCorrection##c.confidence_govt, robust
est store G2

reg YFBehavior YFCorrection##c.confidence_govhealth, robust
est store G3

reg YFBehavior YFCorrection##c.confidence_medicine, robust
est store G4

reg YFBehavior YFCorrection##c.confidence_scientists, robust
est store G5

*Table S24
estout A1 B1 C1 D1 A1b A1c A1d A1e A1f E1 F1 G1, style(tex) replace varwidth(25) collabels("") cells(b(star fmt(%9.2f)) se(par fmt(%9.2f))) stats(N, fmt(%9.2f %9.0f) labels("N")) starlevels(* 0.05 ** 0.01 *** 0.005) 

*Table S25
estout A2 B2 C2 D2 A1b2 A1c2 A1d2 A1e2 A1f2 E2 F2 G2, style(tex) replace varwidth(25) collabels("") cells(b(star fmt(%9.2f)) se(par fmt(%9.2f))) stats(N, fmt(%9.2f %9.0f) labels("N")) starlevels(* 0.05 ** 0.01 *** 0.005) 

*Table S26
estout A3 B3 C3 D3 A1b3 A1c3 A1d3 A1e3 A1f3 E3 F3 G3, style(tex) replace varwidth(25) collabels("") cells(b(star fmt(%9.2f)) se(par fmt(%9.2f))) stats(N, fmt(%9.2f %9.0f) labels("N")) starlevels(* 0.05 ** 0.01 *** 0.005) 

*Table S27
estout A4 B4 C4 D4 A1b4 A1c4 A1d4 A1e4 A1f4 E4 F4 G4, style(tex) replace varwidth(25) collabels("") cells(b(star fmt(%9.2f)) se(par fmt(%9.2f))) stats(N, fmt(%9.2f %9.0f) labels("N")) starlevels(* 0.05 ** 0.01 *** 0.005) 

*Table S28
estout A5 B5 C5 D5 A1b5 A1c5 A1d5 A1e5 A1f5 E5 F5 G5, style(tex) replace varwidth(25) collabels("") cells(b(star fmt(%9.2f)) se(par fmt(%9.2f))) stats(N, fmt(%9.2f %9.0f) labels("N")) starlevels(* 0.05 ** 0.01 *** 0.005) 

reg EPE ZikaCorrection if zika==1, robust 
est store A

*Table S21
estout A, style(tex) replace varwidth(25) collabels("") cells(b(star fmt(%9.2f)) se(par fmt(%9.2f))) stats(N, fmt(%9.2f %9.0f) labels("N")) starlevels(* 0.05 ** 0.01 *** 0.005) 

reg EPE YFCorrection if yellow_fever==1, robust 
est store B

*Table S29
estout B, style(tex) replace varwidth(25) collabels("") cells(b(star fmt(%9.2f)) se(par fmt(%9.2f))) stats(N, fmt(%9.2f %9.0f) labels("N")) starlevels(* 0.05 ** 0.01 *** 0.005) 

reg zika_bite_belief ZikaCorrection##ib4.response_quartile, robust
est store B

reg zika_neuro_belief ZikaCorrection##ib4.response_quartile, robust
est store A

reg zika_sex_belief ZikaCorrection##ib4.response_quartile, robust
est store C

reg zika_handshake_belief ZikaCorrection##ib4.response_quartile, robust
est store D

reg zika_immune_belief ZikaCorrection##ib4.response_quartile, robust
est store E

reg zika_pregnancy_belief ZikaCorrection##ib4.response_quartile, robust
est store F

reg mom_zika_belief ZikaCorrection##ib4.response_quartile, robust
est store G

reg zika_blood_belief ZikaCorrection##ib4.response_quartile, robust
est store H

reg zika_paralysis_belief ZikaCorrection##ib4.response_quartile, robust
est store I

*Table S30b
estout A B C D E F G H I, style(tex) replace varwidth(25) collabels("") cells(b(star fmt(%9.2f)) se(par fmt(%9.2f))) stats(N, fmt(%9.2f %9.0f) labels("N")) starlevels(* 0.05 ** 0.01 *** 0.005) 

reg zika_bite_belief ZikaCorrection##ib4.response_quartile_exp, robust
est store B

reg zika_neuro_belief ZikaCorrection##ib4.response_quartile_exp, robust
est store A

reg zika_sex_belief ZikaCorrection##ib4.response_quartile_exp, robust
est store C

reg zika_handshake_belief ZikaCorrection##ib4.response_quartile_exp, robust
est store D

reg zika_immune_belief ZikaCorrection##ib4.response_quartile_exp, robust
est store E

reg zika_pregnancy_belief ZikaCorrection##ib4.response_quartile_exp, robust
est store F

reg mom_zika_belief ZikaCorrection##ib4.response_quartile_exp, robust
est store G

reg zika_blood_belief ZikaCorrection##ib4.response_quartile_exp, robust
est store H

reg zika_paralysis_belief ZikaCorrection##ib4.response_quartile_exp, robust
est store I

*Table S31b
estout A B C D E F G H I, style(tex) replace varwidth(25) collabels("") cells(b(star fmt(%9.2f)) se(par fmt(%9.2f))) stats(N, fmt(%9.2f %9.0f) labels("N")) starlevels(* 0.05 ** 0.01 *** 0.005) 

